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Summary. — The multiple photon production in Compton scattering at 
high energies is investigated from the Feynman-Dyson formalism of 
quantum electrodynamics. It is observed that the multiple Compton 
scattering cross-section is largest when the photons emitted are all soft 
except one and differs from the single Compton scattering by a multi- 
plicative factor, which increases with increasing energy of the incident 
photon. The multiple cross-section for bremsstrahlung at high energies 
is then obtained by the Williams-Weizsdcker method. The emission of 
about four hard photons within an extremely narrow cone is shown to 
be a plausible estimate. 


1. — Introduction. 


The possibility of multiple production of photons during Compton scat- 
tering, bremsstrahlung and electron-positron annihilation has always been 
overlooked, under the impression that the events are extremely rare owing 
to the feeble interaction in electromagnetic processes. In fact it is obvious 
in the light of the orthodox perturbation theory that probability for the pro- 
duction of » photons will always be smaller by a factor «*1, where 
a = e2JAnhe, than the probability for a similar process involving production 
of a single photon. However, Gupra’s recent calculations on multiple pro- 
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duction of photons in electron-positron annihilation (1) and in brems- 
strahlung (2) showed that at extremely high energies, of the order of 10! eV, 
the effect of the small factor of fine structure constant « was mostly com- 
pensated by another large logarithmic factor which depended on the incident 
energy, and as a result, the cross-section indicated a reasonable probability 
of production of a finite number of photons. These results encouraged us to 
investigate the problem of multiple Compton Scattering, i.e. the scattering 
of a single incident photon into n photons. Calculations are carried out on 
the basis of Feynman-Dyson formalism of quantum electrodynamies (4). It 
is shown that the cross-section is maximum when for a given incident hard 
photon, all the emitted photons are soft, except one, which is hard. The total 
cross-section differs from the single Compton scattering cross-section only by 
a multiplicative factor which depends on the energy of the incident photon 
and vanishes when all the photons are emitted in exactly the forward direction. 

The multiple Compton scattering cross-section is then used to obtain the 
cross-section for multiple bremsstrahlung by the Williams-Weizsàcker method. 
It is shown that emission of up to four photons within an extremely narrow 
cone cannot be ruled out, at energies of the order of 10! eV. Although mul- 
tiple bremsstrahlung cannot be considered an explanation of the Schein 
shower (5), it will, however, be interesting to look for an actual event in which 
a small number of photons, as for instance suggested by our calculations, are 
observed. Discovery of such an event would be a valuable test for the essential 
features of quantum electrodynamics. 


2. — Multiple Compton scattering. 


In the present work we shall throughout use the natural system of units 
h—=c—1 and a metric tensor having non-vanishing components g;1 = 9» = 
= 33 = — Jou = 1. We first consider the multiple scattering of a photon & 
into photon k’, k,, ka, ..., kn, by an electron initially at rest (p= 0, p, =m) 
but finally having a 4-momentum p' (= p’, ps). From the energy-momentum 
conservation relations it can be seen that, for à given incident hard photon, 
only the following set of photon emissions can take place: a) all emitted 
photons are soft, b) one photon (k’) is hard and the rest of the photons are 
soft, and c) more than one emitted photon is hard. Following JOSEPH (°) 


(1) S. N. Gupta: Phys. Rev., 98, 1502 (1955). 

(2) S. N. Gupta: Phys. Rev., 99, 1015 (1955). 

(8) R. P. Feynman: Phys. Rev., 76, 769 (1949). 

(4) F. J. Dyson: Phys. Rev., 75, 486, 1736 (1949). 

(5) M. ScHEIN, D. H. HASKIN and R. G. Grasser: Phys. Rev., 95, 855 (1954). 
(6) J. JosePH: Phys. Rev., 103, 481 (1956). 
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it can be easily shown that the main contribution to the cross-section for the 
production of » photons arises: A) when one of the emitted photons k’ is 
hard, and the rest of the photons k,, k., ..., k, are soft, 
and B) in the Feynman diagrams none of the soft pho- 
tons is emitted by internal emission between the inci- 
dent photon k and the photon k’. It is this case in 
which we shall be interested in the following investiga- 
tion. The Feynman diagram for this case is given beside. 

The matrix element for this graph is given by: 


qn+1(— 4e)"t2 m2 4 È 
(1) Mez 9 ari | ; | OP; — Pr): Fig. 1. — Feynman 
Er) PoPo diagram for multiple 
photon emission by 


v * (p' (20)? P' I k, Mm) Co.) (p' pk + k— m) 1... Compton scattering. 
-€r (p' +k, 4 | k,— my! e' I k k = iy km 
eee (20, )8 P f Ta MACON ip (20) P n n—1 *** r+. 
SR a eT a 
moe Ù à ip si se 2 (20,4) ( kn my kr+s m) oo 
e 
kn pre 
(p m) Co.) u(p) 
where we have put: 
[p=ny, k=k-y, 


(2) 
RIP, PARERI Rai 


Interchanging % and k’ we have a graph G’, the matrix element of which 
is M'. We thus readily find that for |k,|,..., |k,|<|p| and |p'|: 


in(— ie)" 1 DR 
(3) ie Cr} (20, RTE È 
(p' +1) «+. (p''&) i (D 6:41) +. (DP *n) mag 
“(pl hy). (Di + LE) (ph t + ka) (Dn) pe 220 
where, | 
i(— ie)? ( m2 \} 1 : 
A Gn) fa Cw-20) . 


-u(p')[e'(p + k—m)e + e(p + k' — m)te'Ju(p) . 


(4) corresponds to matrix element for the emission of a single photon. We 
are now to sum up all possible matrix elements which arise out of the per- 
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mutations of the soft photons. To do this: A) we first consider the positions 
of photons k and k' fixed on the graph so that r of the soft photons are on 
one side and n —r on the other (a typical example of this is the graph @). 
Then corresponding to this fixed « partition r» we sum over all possible per- 
mutations (p) amongst the soft photons which consist of i) p, permutations 
of soft photons on the p-side; ii) p, permutations of soft photons on the 
p'-side, and iii) interpermutations p,,, between the soft photons of p and 
p'-sides. B) we next sum over all possible partitions 7. This is denoted 


by >. 


Thus we obtain the matrix element for the production of (n+1) photons: 


er 1 
Ru Mot fonts Cam aaah 
(p' -€;) ... (p' -,) (P °C r41) «++ (Pen) 
4 — 1)n- : M,64 i— 
2 ( > (p' +hy)...(p" “hy + ... + Kr) (D Kita t +) (pba) ha Deo 
where, 


2=222. 


(2) Pp Mp’ Tint 
We now use the identity due to GUPTA (1) 


5 1 { 
=, Yul + Va) ce (Yr ++ Ya) Vie Yn 


(6) 


where the summation is taken over all the permutations of the variables, 
and get: 

e” 1 o 
(nr) (2a, ... 20n)} 


(7) Muti Ta 


pin © [PP 6) (Pera) (Dn) A 
ZL Ghia 


(Pi — Di) - 


We next apply the identity due to JOSEPH (°), 


n 


(8) YE) a. dd I] ea)" 
(7) Pint s=1 
and obtain finally: 
eZ TE 1 p'-e, Del mig 
(9) Mai ne Ha (2m,)* ia “ks p Fi] Mò (Di Pi) , 
(10) = Ra dpi — pr). 
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The total cross-section for emission of (n+1) photons is now easily cal- 
culated from MoLLER (7) invariant formula, given by: 


(11) Ont, = (27)? Sar asp! | ask’ | die | dk, > | Rats |?64(p: — pr) - 
pol., spin 

We shall be interested in calculating the differential cross-section for the 
emission of a hard photon k’ into the solid angle dQ’ in the direction 0’, g’ 
and soft photons k,, s= 1,..., n, the photon k, being emitted into the solid 
angle dQ, in the direction 0,, gs. The angles refer to polar co-ordinate system 
with axis in the direction of the incident photon. After integration over 
d*p'dw', which simply removes d-function, replacing p' everywhere in the 
integrand by k—k’—k, ...—k, & k—k', and produces a factor dw'/(d(p,+o’)), 
we obtain: i 
(12) dons, = Cay IT I,w2da, sin 0.00 dp ory 


ihe 


dQ’ > |M}, 


where, 


x \2/ m? ali 
a) X|Mp= (=) (| aay 


and 


p''e pes 2 
pl, p'k 


e? Ai 


(14) DL Ga) Co) È 


The factor n! in the denominator takes into account the fact that the soft 
photons are indistinguishable. Using the polarization summation formula: 


-k\2 
(15) | 2Iprel = (23) 


(00) 
we obtain for p= 0, 


Se iat ON ps i =) 
(16) d = (2a)? al -k.)? = (Be J ’ 


where p’, as mentioned before, stands for k —k’. 
Now from the conservation of momentum 


(17) p -k, = (k —k’)-k, a 


= wo, 089, — w'w, (cos 0, cos 6’ + sin 0, sin 0’ cos (p'— pi)) . 


(7) ©. MoLLeRr: Kgl. Danske Videnskab Selskab, 23; No. 1 (1945). 
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where y’ and g; are the azimuthal angles of the vectors k' and k,. Thus for 
0'<1 we obtain, using conservation of energy: 


(18) p'-ks = pyos(1 — B cos 8.) , 
where, ssa 
i pm _ |p’ 

19 a 7 = te 
al P Do Po 


The expression J, reduces to: 


e? ali b? sin? 0, 
(20) L= af o) cid —f cos 8,)? : 


By integrating over the soft photon variables in (12), we finally obtain the 
differential cross-section for multiple Compton scattering: 


x pe do, sin? 0,40, 
27 ws] (1 — B cos 6,)? 
29 = 1 % Lo tto In & fr Yd I 
(22) Palas TE 7 || Clos 0) du, 


where o, dw’ is the differential cross-section for single Compton scattering 
and is given by: 


dio! 2 ra 
(23) o,do' = arm = È 12 (= n) i 2m | -3) ) 
O) o oO 


@ (60) 


n 


af 
(21) On+1(@, ©) de! = n! TI 
+ 521 


$= 


0,(@, w') do 


and, e and À are the upper and lower limits for the energy of soft photons. 

We see from (21) that the maximum contribution to the cross-section re- 
sults from photons emitted at an angle 0,= cos-!f. Since the factor 
(1/6) In ((1+B)/(1—B)) — 2 in (22) increases as f increases, the cross-section is 
largest as B >1. This corresponds to p>m or w>w'. (22) can then be 
written in the form: 


lle 2p) € 
(24) Onis do! = I, É iù (A) ni (5 


n 
o, do’. 


It is to be emphasized that the multiple Compton scattering is most probable, 
and given by (24), only when all the emitted photons except w’ are soft, 7.¢., 
O DM, Dis Dos On EM, and that @>m, p, >Mm. The multiple Compton 
scattering cross-section differs from the single Compton scattering cross-section 
by a multiplicative factor, and becomes increasingly more important with 
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the increase of incident photon energy. It should be noted that the multiple 
Compton scattering vanishes when the photons are emitted in exactly the for- 
ward direction. A similar result for double scattering has already been dis- 
cussed by MANDL and SKYRME (8). 


3. — Multiple bremsstrahlung. 


We now apply Williams-Weizsacker method (*!) of virtual quanta to 
obtain the multiple photon emission in bremsstrahlung by an electron, moving 
with a high energy £ and velocity v in the field of a nucleus of charge Ze. Let 
us now consider the process in a Lorentz frame of reference in which the 
electron is at rest and the nucleus is moving with a velocity v towards it. 
Now for ve the field of the nucleus is largely transverse, and by Fourier 
analysis it may very well be represented by light quanta travelling parallel 
to the direction of the motion of the nucleus. These virtual quanta suffer 
Compton scattering by the electron at rest and the scattered quanta, on trans- 
forming back to the original system in which the nucleus is at rest and 
electron is moving past by it, appear as the emitted radiation in multiple 
bremsstrahlung. 

The Williams-Weizsäcker spectrum for virtual quanta lying between @ 
and w+dw is given by: 

2 E\ do 
(25) q(@) do = = aZ* In (=) pre 


where x is a constant of the order of unity. The differential cross-section for 
multiple bremsstrahlung is then given by: 


(26) o (odo! = do’ to) dos. 0, DIG 


where 051(0, w') do’ is the cross-section for multiple Compton scattering given 
by (21). Transforming back to the laboratory system and integrating we 
obtain the total cross-section as: 


(27) De. É In Fa In (5) 


N. |IT 


(8) F. ManpL and T. H. SkyRME: Proc. Roy. Soc., A 215, 497 (1952). 
(2) E. J. Wixriams: Kgl. Danske Videnskab Selskab, 13, No. 4 (1935). 
(19) C. F. v. WEIZSACKER: Zeits. f. Phys., 88, 612 (1934). 
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where b is a constant of the order of unity, and o? is the corresponding Bethe- 
Heitler bremsstrahlung formula given by: 


«A 


We now note the following relations between 6 and w and their Lorentz trans- 
forms ò,, 0, in the laboratory system; 6 is the direction of the recoil electron 
with respect to the incident photon. 


(28) oF = 47°07? 


For p,>M: 
| Ri ne SUR 
ted = tg = ee OR) 
(29) and 
Lara Be 
Re ae i =F 


Since in Compton scattering the electron cannot be scattered backwards, it 
is clear that the deflected electron goes otf at an angle 6, > m/E with the 
incident electron direction. The energetic photons are emitted in a narrow 
cone around the direction of the incident electron, while the softer photons 
are emitted at larger angles. It is therefore only those photons which are 
emitted in a narrow cone which will be susceptible to observations. We will 
therefore be interested in finding out the cross-section for multiple brems- 
strahlung of hard photons lying within a narrow cone between angles 0 and 6, 
around the direction of the incident electron. We therefore carry out the 
integration in (21) over angles 0, lying between 2m/9, E and x instead of between 
0 and x, as was done earlier, and obtain: 


x ai | 147 | 
30 Pal) — 
EEE 
Il 1 e\|" 
ei paz 1 o 
DE) ira n(3) 4 
where, 
Mm 2 
(81) n =) 
This reduces to 
SI 1 [2 Eô, é\|" 3 
(32) Ce = El a | In i] Gia for ae 
and, 5 
È 1 [2x (E60, svi" 
(33) thes = fe in (7) mm (5) di, for DORE 
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It is to be observed that the cross-section for multiple Bremsstrahlung is 
maximum at 


Po — M 
(34) cos 0, = B = si 
giving, 
2m\} 
35 0, RI (=) . 
(35) nl 


When transformed to the laboratory system by the relation (29) the maximum 
occurs at 0, & 0,. In order that the maximum may be within the range of 
angular integration from 0 to 6, we must have 0,> 6, and the cross-section 
which is of interest for multiple Bremsstrahlung is given by (32). 

Now we calculate the order of magnitude of multiple bremsstrahlung as 
compared with the single Bremsstrahlung. Taking E —10! eV, 6, ~ 10-4 radians 
e/A — 108, we find (2x/7) In (Ed, /m) In (e/A) & 1, and the cross-section becomes 


i 
(36) One® "i Ge 


We therefore see that the emission of about four hard photons within a 
narrow cone cannot be ruled out. Experimental realization of such events 
will give a valuable check on the validity of quantum electrodynamics. 


Two of us (V.S.M. and J.D.) are thankful to the University Grants Com- 
mission for the award of fellowships which made this investigation possible. 


RIASSUNTO (*) 


Partendo dal formalismo di Feynman-Dyson dell’elettrodinamica quantistica si 
studia la produzione multipla dei fotoni nello scattering Compton alle alte energie. 
Si osserva che la sezione d’urto per lo scattering Compton multiplo è massima quando 
i fotoni emessi sono tutti molli eccetto uno e differisce da quella dello scattering Compton 
semplice per un fattore moltiplicativo che cresce coll’aumentare dell’energia del fotone 
incidente. La sezione d’urto multipla per la bremsstrahlung alle alte energie si ottiene 
poi col metodo di Williams-Weizsäcker. Si mostra che l’emissione di circa quattro 
fotoni duri entro un cono estremamente ristretto risponde a una previsione plausibile. 


(*) Traduzione a cura della Redazione. 
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On the Determination of the Parity of Hyperons and K-Mesons. 


S. M. BILENKY and R. M. RYNDIN 


Joint Institute for Nuclear Research, Laboratory of Theoretical Physics - Dubna 


(ricevuto il 29 Novembre 1958) 


Summary. — Possible methods for determining the relative parity of 
hyperons and K-mesons in the reactions with polarized proton beams were 
considered. 


The determination of the intrinsic parity of K-mesons and hyperons is 
one of the most important problems of the physics of elementary particles. 
Since the parity is not conserved in all weak interactions including the decays 
of strange particles the intrinsic parity of hyperons and K-mesons may be 
determined only in the study of strong interactions responsible for their pro- 
duction and mutual transformations. 

Recently some ways for the determination of the parity of strange particles 
.based upon the study of the reactions involving hyper-nuclei and of those 
involving strange particles near the threshold (**) have been suggested. 

Other possibilities are pointed out in this note. 


1 


) R. H. DALITZ: Rep. on Progr. in Phys., vol. 20 (London, 1957). 
2) CHou Kuane-Cuao and M. I. Srroxoy: Nucl. Phys., 6, 10 (1958). 
3) CHou KuaxG-CH40: Journ. Exp. Theor. Phys., 34, 1027 (1958). 

) L. B. Oxun and I. J. POMERANCUK: Journ. Exp. Theor. Phys., 34, 997 (1958) 
5) G. Costa and B. T. FeLD: Phys. Rev., 109, 606 (1958). 

(6) L. B. Oxun, I. J. PomERANCUK and I. M. SmuskEvié: Journ. Exp. Theor. 


Phys., 34, 1246 (1958). 
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1. — Consider the reaction (*) 


(1) p+ ‘He > ‘He4+- Y4 K 
with a polarized proton beam. The matrix of reaction (1) is of the form 
(2) M = a—+o:b, 


where the quantities a and b are functions of the momenta of the initial 
and final states. The spins of hyperon and K-meson are assumed to be } 
and 0 respectively. Using usual methods (7) we obtain the following expression 
for the differential cross-section of reaction (1) involving the polarized proton 
beam 


9 * ey 
(3) o=o(1 pee b + ib “i 


aa” + b-b* 


Here o, is the cross-section of process (1) for unpolarized incident beam, 
whereas P, is the initial proton polarization. On the other hand, hyperon po- 
larization in reaction (1) with the unpolarized protons is 


2Rea*b + ib x b* 
aa* + b - b* 


(4) PS 


The product of the intrinsic parities of proton, hyperon and K-meson may 
take two values I, J,J,=-+1. In the first case the matrix M is a scalar, 
i.e., the quantity b is an axial vector and may be presented in the form 


(5) b = b, Kx K’+ (K- K’x K")(b, K+ b,K’), 


where K, K' and K” are the momenta of a proton, hyperon and K-meson (**) 
respectively. The quantities a and d; are arbitrary functions of scalar products 
constructed from the vectors K, K’ and K”. 

In the second case M is a pseudoscalar, i.e. the quantity b is a vector 


(6) b = (K-K'x K’)b, Kx K'+b,K+b,K'. 
(‘) Note that in (1) helium may be replaced by any nucleus with zero spin, e.g. °C. 


(7) L. WoLFENSTEIN and J. ASHKIN: Phys. Rev., 85, 947 (1952). 
(**) Here and further on the momenta K and K' are not parallel to each other. 
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The quantity a is a pseudoscalar 

(7) a = (K-K’'x K")a’. 

a’ and b, are arbitrary scalar functions. De. 


If the momenta of proton, K-meson and hyperon lie in the same plane 
then the triple scalar product (K-:K'x K”) vanishes and if 1,1,1,=1 


(5’) b=0),Kx K', 
Whale if LT Te = 1 

(6’) b = b, K+), K’ 
and 

(7°) ai 


Then, as can be easily seen from (3) and (4), the expression for the differential 
cross-section takes the form 


(8) o=0o(1 + PP), 


where the sign is (+) if the complete intrinsic parity does not. change, and 
(—) otherwise. 

Polarized proton beams are obtained in the scattering of unpolarized beams 
by nuclei and their polarization is orthogonal to the momentum of the scat- 
tered protons. Choosing the z-axis in a direction of vector K, and the y-axis 
in that of P,, we write (8) as follows 


(9) o=o0(14+P.P cos y’). 
Hence, the azimuthal asymmetry of the hyperon distribution is 
(10) : gie + P,P’ 
Thus, the measurement of the azimuthal asymmetry in reaction (1) with 
polarized protons makes it possible to define the sign of I,IyIx uniquely, 


if the sign of the polarization P arising as a result of the reaction with the 
unpolarized protons (*) is known. The sign of P may be determined by 


(*) The magnitude and sign of the polarization P, are supposed to be known from 
previous experiments on double proton scattering. 
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studying the pion asymmetry from the decay of polarized hyperons. As it 
follows from the previous arguments in the determination of the azimuthal 
asymmetry it is necessary to choose such events in which all the particles lie in 
the same plane. One may perform the integration over the angle between 
the directions of the momenta of a proton and K-meson. 

Another reaction of the considered type is 


(11) p+‘He >‘H,+p+K*+. 


if the spin of a hypernucleus *H, is zero. However, reaction (1) considered 
above does not require such an assumption. Moreover, its study would enable 
to determine the relative parity of different hyperons. 


2. — Consider the reaction (3) 


(12) p+*He >*He, + Kr. 


Suppose the spin of the hypernucleus *He, to be 3. Since in this case only 
two vectors (the relative momenta of the initial and final states) are avail- 
able then, as is seen from (2)-(4), the cross-section of the reaction (12) 
with the polarized proton beam is of the form (8). Therefore the measu- 
rement of the azimuthal asymmetry of K-mesons makes it possible to deter- 
mine the relative parity (AK) if the sign of polarization of the hyper-nucleus 
*He, in reaction (12) with the unpolarized protons is known. These conclu- 
sions refer also to the reactions with a polarized antiproton beam. For 
example, to the reaction 


(13) p+‘He > *H, + Ke 
if the spin of *H, is 3. 

In conelusion let us point out that reactions (1) and (11) when all the 
particles lie in the same plane as well as reactions (12) and (13) are processes 


of the type 


(14) Pot. 


For such processes the differential cross-section always has the form (8), 
that makes it possible to use them for the unique determination of the internal 
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parity. Here there are two possibilities: 
1) Polarized target with spin 3 and a beam of spinless particles (°°). 


2) Polarized beam of particles with spin 3 and a spinless target. 


The latter possibility has been considered in thi note. 


The authors are grateful to CHOU KUANG-CHAO, B. PONTECORVO and YU. 
SCERBAKOF for interesting discussions of the results obtained. 


(8) S. M. Birenxy and R. M. RyNDIN: Journ. Exp. Theor. Phys., 35, 827 
(1958). S. M. Birenxy: Journ. Exp. Theor. Phys., 36, 291 (1959). 
(9) S. M. Brrenxy: Nuovo Cimento, 10, 1049 (1958). 


RIASSUNTO (*) 


Si esaminano i possibili metodi per la determinazione della parità relativa degli 
iperoni e mesoni K nella reazione con fasci di protoni polarizzati. 1 


(*) Traduzione a cura della Redazione. 
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Le caratteristiche della macchina definitiva del C.S.C.E. di Pisa 
dal punto di vista logico-matematico (*). 


A. (CARACCIOLO (*) e L. GUERRI (*) 


Centro Studi Calcolatrici Elettroniche dell’Università - Pisa 


(ricevuto il 12 Dicembre 1958) 


Riassunto. — Si descrivono le caratteristiche generali della calcolatrice 
in costruzione presso il C.S.C.E. con particolare riferimento alla modifica 
automatica degli indirizzi mediante celle parametriche. 


1. — Introduzione. 


La calcolatrice elettronica in costruzione presso il C.S.C.E. (Centro Studi 
calcolatrici elettroniche) di Pisa, ha le seguenti caratteristiche: è in « parallelo » 
con parole di lunghezza fissa di 36 bits. Una parola rappresenta un numero 
o una istruzione. I numeri sono rappresentati, in base 2, in virgola fissa o in 
virgola mobile. Per la virgola mobile 28 bits sono per la parte frazionaria, 8 
per l'esponente. La memoria principale a nuclei magnetici, sarà di 4096 parole, 
ma è prevista la possibilità di ampliare la memoria fino a 32768 parole senza 
apportare modifiche alla calcolatrice. La memoria ausiliaria sarà costituita da 
un Tamburo Magnetico di 16384 parole e da alcune unità a nastro magnetico. 
Il numero dei nastri magnetici collegati con la calcolatrice quando questa 
entrerà in funzione non è stato ancora fissato ma, come nel caso della memoria 
principale è prevista la possibilità di aggiungere successivamente varie unità 
a nastro magnetico. 


(*) Presentato al XLIV Congresso della Società Italiana di Fisica, Palermo, 6-11 


Novembre 1958. 4 
(**) Distaccato dalla Sezione di Pisa dell’I.N.F.N. presso il C.S.C.E. 
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2. — Istruzioni. 


La calcolatrice del Centro ha una lista di oltre 100 istruzioni. Ogni istru- 
zione è rappresentata con una parola. Le istruzioni sono del tipo ad un indi- 
rizzo con modifica e si distinguono in ordinarie e speciali. Questa distinzione 
è fatta in base al diverso uso delle celle parametriche (si veda ‘più avanti per 
la definizione di cella parametrica). Le istruzioni ordinarie sono modificabili 
fino a due volte, quelle speciali una volta. Una caratteristica della calcolatrice 
del Centro è la soluzione adottata per la modifica delle istruzioni che è fatta 
utilizzando celle ordinarie della memoria invece di registri speciali. 


3. — Celle parametriche. 


Oltre ai registri aritmetici di cui si fa cenno in altra parte della nota, la cal- 
colatrice ha due registri, indicati con H° e H!, di 15 bits ciascuno. Per mezzo 
di una istruzione si può registrare in H° (in H!) un numero intero 4° (41). Le 
celle di indirizzo: 


RL ROSES NN RAIN e I 


formano due gruppi di celle che si indicano con {H°} e {H+}. Le celle di questi 
gruppi si chiamano celle parametriche. L'indice di gruppo che vale 0 o 1 speci- 
fica uno dei due gruppi. I numeri 1, 2 ... 31, si dicono indirizzi relativi delle 
celle parametriche. 

h° (h1) si dice origine di {H°} ({H*}). 

Fissati 4° e h! una cella parametrica è completamente individuata per mezzo 
dell’indice di gruppo e dell’indirizzo relativo. Cambiare il contenuto del registro 
H significa considerare da quel momeno in poi e fino a nuovo cambiamento, 
come gruppo {H°} le 31 celle consecutive a quella che ha come indirizzo il nuovo 
contenuto di H°, Analogamente per H1, 

Le celle parametriche svolgono lo stesso ruolo di registri di modifica nel 
modo che: viene descritto più avanti. Le celle parametriche di indirizzo rela- 
tivo 0, cioè indirizzo assoluto 4° e h1 non fanno parte del gruppo. Dal punto di 
vista del programmatore è come se tali celle contenessero 0 e il loro contenuto 
non potesse venir cambiato. 

È importante osservare, che le celle parametriche, comprese quelle di indi- 
rizzo relativo 0, sono nello stesso tempo disponibili come ordinarie celle di me- 
moria senza restrizioni speciali. 

In una istruzione le celle parametriche sono usate nel modo seguente: siano 
Sos S1 --. 835 1 36 bits dell’istruzione. Essi sono suddivisi: 


(So, 81) (Sz +++ 88); (Sg +++ S14)(S15 +++ 820) (821 «++ 835) + 
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So: Sì Sono usati per arresti controllati dal quadro di comando. 

8: .. Sg danno il simbolo di operazione. 

So © S10 -- Su danno rispettivamente l’indice di gruppo e un indirizzo rela- 
tivo che definiscono una cella parametrica P’. 

815 © Sig + Soo danno rispettivamente un indice di gruppo e un indirizzo rela- 
tivo che definiscono una cella parametrica P’. 

So, «+. 835 definiscono un numero intero compreso fra 0 e 32 767. 

Con P' e P” si indicano quindi le celle parametriche definite da (ss .. 814) e 
(815 ++ S20)- 

Indichiamo con s il numero intero definito da (8; ... s3;)e con p' (p’) il conte- 
nuto di P’ (P”) con l’intesa che p'—0 (p"=0) se Vindirizzo relativo P’ (P") è 0. 

Nelle istruzioni ordinarie, l’istruzione è eseguita con Vindirizzo s+p'+p’. 

Così se l’istruzione ordinaria è modificata una volta sola è indifferente con- 
siderare la cella impiegata nella modifica come P' o P”. 

Allo stesso modo se le modifiche sono due è indifferente l’ordine in eui si 
considerano le celle parametriche usate. 

Invece, per quanto riguarda le istruzioni speciali P’ e P” hanno ruolo di- 
Verso. 

P' è usato nello stesso modo descritto per le istruzioni ordinarie. 

P" invece può essere usato in diversi modi a seconda dell’istruzione spe- 
ciale che si considera. 

Indichiamo con A(s+p’) e yA(s+p') rispettivamente la cella di indirizzo 
s + p' e il contenuto di tale cella. Vi sono per esempio le seguenti istruzioni 
speciali: 


YT salti condizionati o incondizionati rientrano tra le istruzioni speciali. 
L'indirizzo dell’istruzione a cui si salta è dato da s+4+p' e in P” viene inviato 
l’indirizzo presente nel contatore di istruzione al momento dell’esecuzione del 
salto. Questo fatto facilita il rientro dai sottoprogrammi nel programma prin- 
cipale. 

In definitiva nella calcolatrice del Centro sono disponibili due gruppi di 31 
celle l’uno che funzionano come registri di modifica. Questi gruppi si possono 
sostituire con altri con un’apposita istruzione che definisce l’origine dei nuovi 
gruppi. 

Un inconveniente di questo metodo è che ogni modifica richiede una lettura 
della memoria e quindi circa 10 us supplementari. Non si ha perdita di tempo 
se l’indirizzo relativo è zero. 


8 - Il Nuovo Cimento. 
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4. — Registri aritmetici. 


I registri aritmetici a cui fanno riferimento le istruzioni sono 2 indicati con 
A e B. Ognuno di questi registri è di 36 bits. B ha un ruolo ausiliario nei con- 
fronti di A ed è utilizzato in particolare per la moltiplicazione, la divisione e 
i calcoli di precisione multipla. 

Un’altra caratteristica della Calcolatrice è la possibilità di eseguire alcune 
operazioni tra unità aritmetica e memoria con registrazione del risultato in A 
oppure in memoria. Così, indicando con a il contenuto di A, e con 2 il contenuto 
di una cella di memoria Z, sono possibili le seguenti operazioni : 


g > À 
—g> À 
2 a VA 
— 4 —+ À 


ASA 


2 
| 
2 


con Z immutato 
con Z immutato 
con Z immutato 
con Z immutato 


con Z immutato 


a—Z con A immutato, 
—a—->Z con A immutato , 
a+z—Z con A immutato, 
è—a—Z con A immutato , 


a—&->Z con A immutato, 


nelle operazioni aritmetiche il traboccamento è sempre segnalato su apposito 


indicatore. 


Per quanto riguarda i salti oltre il salto incondizionato, vi sono i seguenti: 


salta 
salta 
salta 


salta 


sea>0, 
se a<0, 
se a=0, 


se a 0, 


a salti analoghi su B. Inoltre vi sono salti su traboccamento. Vi sono pure 
istruzioni di confronto e istruzioni di salto su confronto. Tra queste istruzioni 
di confronto vi sono per esempio: 


confronta a e 2, 


confronta ze 0, 


confronta i-mo bit di 2 con 0. 


Nel confronto a e + rimangono immutati. 
I salti su confronto sono: 


— salta se nell’ultimo confronto i due numeri erano =; 


— salta se nell’ultimo confronto i due numeri erano 
— salta se nell’ultimo confronto il primo era > del secondo; 
— salta se nell’ultimo confronto il primo numero era < del secondo; 
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La moltiplicazione in virgola fissa ha tre varianti: 

— moltiplicazione cumulativa con risultato di 72 bits in A e B; 
— moltiplicazione esatta con risultato di 72 bits in A e B; 
— moltiplicazione arrotondabile con risultato di 36 bits in A. 

La divisione in virgola fissa ha tre varianti: 

— divisione lunga esatta: il dividendo di 72 bits ha i primi 36 bits in A e i 

rimanenti in B. Il quoziente si forma in A il resto in B, 

— divisione lunga arrotondabile: il dividendo è in A e B, il quoziente si forma 
in A; 
— divisione corta arrotondabile: il dividendo di 36 bits e in A il quoziente si 

forma in A. 

Diciamo operazioni con risultato. arrotondabile invece che arrotondato 
perchè con un’apposita istruzione si può controllare l'arrotondamento in tutte 
le successive operazioni con risultato arrotondabile. 

In tal modo uno stesso calcolo può essere ripetuto con o senza arroton- 
damento. 

Questo fatto può fornire qualche indicazione sulla precisione dei risultati. 
Nelle operazioni in virgola mobile l’arrotondamento e la normalizzazione dei 
risultati sono controllabili separatamente nello stesso modo descritto per la 
virgola fissa. 

Riportiamo ora i tempi di esecuzione delle principali istruzioni. Questi 
tempi vanno aumentati di 10 us per la chiamata dell’istruzione stessa; ogni 
modifica dell’istruzione richiede altri 10 us. 

— Trasferimenti tra registri e memoria: 10 us. 

— Somme e sottrazioni in virgola fissa: 10 us. 

— Moltiplicazioni in virgola fissa: 240 us. 

— Divisioni in virgola fissa: 280 us. 

— Somme e sottrazioni in virgola mobile: (80-+4n) us, con n=a numero di 
traslazioni necessarie per allineare la virgola e per normalizzare il risultato. 

— Moltiplicazione in virgola mobile: 210 us. 

— Divisione in virgola mobile: 260 us. 

Questi tempi rappresentano solo una indicazione di massima. 

Sulla base dell’esperienza fatta sulla calcolatrice ridotta già costruita al 
Centro (parole di 18 bits, memoria di 1024 parole, 30 istruzioni senza modi- 


x 


fica) non è esclusa la possibilità di ridurre i tempi dati sopra del 20% circa. 


SUMMARY 


A description is given of the computer being built.at the C.S.C.E. with particular 
reference to the solution adopted for the modification of the instructions by means 
of «parametric cells ». 
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IL NUOVO CIMENTO Vor. XII, N. 2 16 Aprile 1959 


Memorie ed entrata-uscita 
della « macchina definitiva» del C.S.C.E. (*). 


A. CARACCIOLO (**), G. CECCHINI (**), G. B. GERACE (**), 
M. FALLENI (*) e A. SABBADINI (**) 


Centro Studi Calcolatrici Elettroniche dell’ Università - Pisa 


(ricevuto il 28 Dicembre 1958) 


Riassunto. — Si descrive la connessione delle memorie e dell’entrata-uscita 
con il nucleo centrale della macchina in costruzione presso il Centro di 
Pisa. Si illustra il complesso di tali apparecchiature di cui verrà dotata 
inizialmente la macchina e le notevoli possibilità di ampliamenti già 
previsti dal progetto logico generale. Si accenna alla possibilità di colle- 
gamenti diretti con istituti di ricerca interessanti all’uso della macchina. 


1. — Introduzione. 


Il giusto proporzionamento delle capacità e velocità delle memorie di un 
calcolatore elettronico e delle velocità degli organi di entrata e uscita dipende 
da una larga serie di fattori. 

Tra i principali indichiamo: 

— La velocità base di calcolo del nucleo centrale. 

— Il tipo del calcolatore (per usi scientifici o commerciali). 

— Le caratteristiche dei calcoli da eseguire. 

— La quantità dei dati da introdurre in entrata e stampare in uscita. 


(*) Presentato al XLIV Congresso della Società Italiana di Fisica, Palermo, 6-11 
Novembre 1958. 

(**) Distaccato dalla Sezione di Pisa dell’I.N.F.N. presso il C.S.C.E. 

(**) Distaccato dalla S.p.A. Olivetti presso il C.S.C.E. 
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— La struttura logica della macchina e l’organizzazione della programmazione. 


— I costi. 


Nella presente comunicazione, anche se talvolta sarà fatto riferimento al 
problema delle scelte, verrà essenzialmente illustrato il complesso di apparec- 
chiature di cui sarà dotata inizialmente la « macchina definitiva » in costru- 
zione presso il Centro di Pisa e le notevoli possibilità di ampliamento già pre- 
viste dal progetto logico generale. 

Si descriverà la connessione delle memorie e dell’entrata-uscita con il nucleo 
centrale e verrà fatto cenno alla possibilità di collegamenti diretti con istituti 
di ricerca interessati all’uso della macchina. 


2. — Previsioni di progetto. 


Il progetto della « macchina definitiva » prevede: 


1) Una «memoria veloce » a nuclei magnetici che riceve e trasmette le 
informazioni (istruzioni e dati) da e verso il nucleo centrale e le « apparec- 
chiature ausiliarie » del calcolatore. 

2) Una serie di «apparecchiature ausiliarie » che si compone di: a) me- 
morie ausiliarie; b) organi di Entrata e Uscita. 

a) Le «memorie ausiliarie » comprendono un tamburo magnetico ed 
una serie di nastri magnetici (da 1 a 20) che ricevono e trasmettono le informa- 
zioni da e verso la « memoria veloce ». 


b) Gli organi di entrata e uscita comprendono: 


— Uno o più lettori fotoelettrici a 400 caratteri/s in entrata per la lettura su 
banda perforata. © | 

— Una telescrivente in entrata per l’introduzione diretta di dati «numerici. 

— Uno stampatore parallelo di uscita a 250 caratteri/s per la composizione 

e stampa dei risultati. 

— Una telescrivente ricevente di uscita con funzioni diverse tra cui il con- 
trollo di calcoli e programmi con dati significativi. 

— Un perforatore di banda veloce di uscita a 60 caratteri/s per la reintroduzione 
sui lettori fotoelettrici dei dati e programmi perforati. 

— Un quadro di comando manuale per le operazioni di avviamento della mac- 
china ed i necessari controlli. 

I nastri magnetici, oltre che come « memorie ausiliarie », possono essere 
usati come organi veloci di-entrata e di uscita, ammesso di disporre fuori della 
macchina delle opportune apparecchiature commerciali di conversione con i 
normali mezzi perforanti e stampanti. 
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3. — La memoria veloce. 


La «memoria veloce» avrà inizialmente 4096 celle, ma il progetto della 
macchina ne prevede l'ampliamento fino-a 32768, con la semplice aggiunta 0 
sostituzione dei gruppi di memoria. 

Ricorrendo alle tecniche di indirizzamento indiretto connesse all’uso delle 
celle parametriche, è possibile elevare la capacità della memoria a valori prati- 
camente illimitati. 

La tendenza ad aumentare la capacità delle « memorie veloci » è principal- 
mente dovuta alla crescente necessità di risolvere problemi con grande numero 
di dati accessibili a caso (tipico ad esempio l’impiego del metodo di Monte 
Carlo nella risoluzione dei problemi nucleari), alla agevolazione che introduce 
nei moderni metodi di programmazione automatica ed al costo dei nuclei fer- 
ritici che nel giro di pochi anni si è ridotto di molti fattori. 

Attualmente, a parità del numero di celle, una memoria veloce a nuclei 
ferritici viene a costare solo circa il doppio di una memoria lenta a tamburo. 


4. — Il controllo esterno. 


Le « memorie ausiliarie » (tamburo e nastri magnetici) e gli organi di En- 
trata sono forniti di un controllo esterno molto centralizzato e separati dal nu- 
cleo centrale e dalla memoria veloce da un unico Registro Esterno. 

Per l’uscita è utilizzato un registro aritmetico della macchina. La trasmis- 
sione di una parola (36 bit) tra le Apparecchiature Ausiliarie ed il Registro 
Esterno avviene con 6 trasferimenti successivi di un carattere di 6 bit utili. 

Questa soluzione, oltre a permettere di utilizzare istruzioni pressochè iden- 
tiche nell’impiego di alcune Apparecchiature Ausiliarie (ad es, nelle istruzioni 
di trasferimento a blocchi del tamburo magnetico, nastri magnetici e lettori 
fotoelettrici) conduce a poter utilizzare per l’entrata un carattere di 6 bit che 
facilita notevolmente l'introduzione di dati a-numerici per la programmazione 
automatica. 

In realtà, all’interno delle « apparecchiature ausiliarie », ai 6 bit utili viene 
aggiunto un bit per il controllo di parità che nei trasferimenti dal nucleo centrale 
e memoria veloce verso le Apparecchiature Ausiliarie viene aggiunto automa- 
ticamente e nei trasferimenti inversi controlla il complesso delle operazioni 
serittura-lettura. 

L'introduzione di un bit di controllo solo nelle « Apparecchiature Ausiliarie » 
è giustificata dalla loro minore sicurezza per la presenza di organi meccanici 
in movimento. 
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5. — Caratteristiche e istruzioni delle apparecchiature ausiliarie. 


Daremo ora una breve descrizione delle caratteristiche delle Apparecchiature 
Ausiliarie adottate e delle istruzioni che le riguardano. 

Tutte le istruzioni che seguono sfruttano a fondo le possibilità offerte dalle 
celle parametriche e la flessibilità logica del calcolatore. 

Il Tamburo magnetico ha le seguenti caratteristiche: 16384 celle di me- 
moria; 40 us di tempo di trasferimento; 10 ms di tempo medio di accesso per 
un programma non ottimizzato temporalmente. 

Dispone di 5 istruzioni: 


— Una, per l'indirizzo della prima parola sul tamburo da cui iniziare il tra- 
sferimento. 


— Due per il trasferimento nel senso Tamburo magnetico + memoria veloce 
e viceversa di un determinato numero di parole a partire da una determinata 
cella di memoria. 


— Due per il trasferimento Tamburo Magnetico + memoria veloce e viceversa 
di un blocco di parole con in testa la sua lunghezza; e ciò a partire da una 
determinata cella di memoria. 


Il trasferimento Tamburo Magnetico + memoria veloce avviene in memoria 
veloce a cominciare dalla prima parola utile del blocco mentre il numero di 
parole contenuto in un blocco può essere conservato in una cella parametrica. 

Da notare la grande flessibilità del sistema che offre la possibilità di trasfe- 
rimenti a numero di celle ed a blocchi variabili. 

Il Nastro Magnetico (N.M.) ha le seguenti caratteristiche: 250000 celle di 
memoria; 60 us di tempo di trasferimento; 2 ms di tempo d’accesso conside- 
rando il blocco indirizzato già posizionato. 

Dispone di 3 istruzioni: 


— Una per la ricerca del blocco di un nastro determinato. 


La ricerca si svolge mentre il calcolatore seguita ad operare e l'istruzione 
dà contemporaneamente l’indirizzo dell’unità a nastro e quella del blocco. 


— Due istruzioni di trasferimento a blocchi variabili nei due sensi, del tutto 
simili a quelle descritte per il Tamburo Magnetico. 


Il sistema anche qui si mostra particolarmente flessibile confrontato a quello 
di altri calcolatori che difficilmente realizzano contemporaneamente il funzio- 
namento a blocchi variabili e la ricerca indipendente del blocco. 

Per un confronto diretto tra le velocità di impiego del Tamburo Magnetico 
e del Nastro Magnetico osserviamo quanto segue: 
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ammesso un blocco medio di 100 celle e non considerando il blocco in- 

dirizzato già posizionato, abbiamo un tempo di accesso variabile da 2 ms 
a 208. 

Mentre però con un programma ottimizzato spazialmente (quando ad esempio 
il nastro è usato come libreria di sotto programmi) limite superiore si riduce 
notevolmente, con una ottimizzazione spazio-temporale ed in programmi par- 
ticolari la velocità di impiego del nastro può essere vantaggiosamente confron- 
tata con quella del Tamburo Magnetico. 

Per l’Entrata-Uscita abbiamo ancora 4 istruzioni. 


— Una istruzione per l’entrata di un carattere della Telescrivente o da uno 
dei Lettori fotoelettrici. 


— Una istruzione per l’uscita di un carattere verso lo stampatore parallelo, 
il perforatore veloce o la telescrivente. 


In ambedue, l’indirizzo dell’istruzione determina l’organo prescelto in En- 
trata e in Uscita dando la possibilità di aggiungere un numero sufficientemente 
grande di nuove apparecchiature essegnando ad esse nuovi indirizzi. 


— Una istruzione per la lettura a blocchi dal Lettore fotoelettrico per l’intro- 
duzione diretta di programmi senza conversione. 


— Una istruzione di stampa per lo stampatore parallelo. 


Questa istruzione permette la stampa dei dati «-numerici mentre il calco- 
latore seguita ad operare. Ciò viene ottenuto corredando l’organo stampante 
di una memorietta ausiliaria da riempirsi a velocità di macchina con l’istruzione 
di uscita. 

Il Quadro di comando manuale è stato studiato in modo da rendere il più 
possibile semplici e rapide le operazioni manuali dell’operatore. 

Esso è provvisto inoltre di una serie di indicatori e di dispositivi per age- 
volare le operazioni di controllo tecnico della macchina. 


6. — Considerazioni generali. 


Durante l’impiego delle apparecchiature ausiliarie il controllo permette la 
sincronizzazione del calcolatore sulle diverse velocità delle apparecchiature 
ausiliarie. 

Se sopraggiunge un’istruzione riguardante un organo esterno ancora’ impe- 
gnato da una istruzione precedente (organo «non pronto ») il calcolatore cessa 
di operare ed attende che questi divenga « pronto ». Dopo di che automatica- 
mente torna ad operare. 

Non possono essere impegnate contemporaneamente più apparecchiature 
ausiliarie ad eccezione del periodo di stampa per lo stampatore e della ricerca 
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di blocco per il Nastro Magnetico. (In questo caso non possono essere effettuate 
ricerche contemporanee su nastri diversi.) 

È previsto di utilizzare i tempi morti delle Entrate ed Uscite veloci (tempo 
di spostamento di un carattere da una posizione all’altra di lettura per il let- 
tore fotoelettrico e tempo di stampa per lo stampatore parallelo) per il cam- 
biamento di rappresentazione da decimale codificato a binario puro, da ese- 
guirsi su programma all’interno del nucleo centrale. 

Non è esclusa la possibilità di effettuare a parte il cambiamento di rappre- 
sentazione da mezzi perforanti su Nastro Magnetico e viceversa impiegando 
un altro piccolo calcolatore (ad es. la Macchina Ridotta). 

In tal caso, utilizzando come Entrata e Uscita i nastri magnetici del 
calcolatore in progetto, si verrebbe ad operare in Entrata e Uscita ad una ve- 
locità media di circa 70000 caratteri al secondo nettamente confrontabile con 
la velocità di base del calcolatore stesso. 


7. Cenni sul problema dei collegamenti a distanza. 


Per un grande calcolatore destinato a servire Istituti Scientifici diversi e 
lontani sorge il problema di trasmettere rapidamente a distanza dati e pro- 
grammi. 

Mentre per i dati il problema è evidente, nel caso di programmi possono 
darsi due eventualità. 

Una in cui il ricercatore lontano, dovendo risolvere problemi particolarmente 
complessi, debba raggiungere il centro di calcolo per i contatti necessari. 

L’altra in cui i programmi siano affrontabili direttamente dal ricercatore 
con la programmazione automatica. 

È in quest’ultimo caso che sorge la necessità della prova rapida a distanza di 
pezzi di programmi con ricezione altrettanto rapida dei risultati. 

Il problema può essere affrontato con i normali mezzi di comunicazione ma 
esiste anche la possibilità di risolverlo con l’impiego dei mezzi di telecomuni- 

cazione. 

Negli Stati Uniti si ha già una buona esperienza su problemi del genere, 
ed essi vengono soddisfacentemente risolti con trasmissione diretta su tele- 
scrivente a 5 canali, considerata la più esente da errori. 

Anche lo standard Italiano per le linee e le apparecchiature funziona su 5 
canali non controllati mentre l’ingresso del nostro calcolatore è a 7 canali, di 
cui uno di controllo. | 

Il problema può essere facilmente risolto fuori della macchina con una appa- 
recchiatura di conversione di codice che trasformi la banda ricevente da 5 a 
7 canali e viceversa; oppure aggiungendo alla macchina un’entrata-uscita a 
cinque canali corredata da opportuno programma di conversione. 
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Il controllo sulla trasmissione di linea viene sostanzialmente risolto con la 
trasmissione doppia della banda ed il controllo dell’equivalenza all’arrivo. 

Per l’attrezzatura dei Centri lontani è sufficiente una telescrivente-rice-tra- 
smittente munita di ‘perforatrice ed un trasmettitore automatico del costo 
complessivo di 500000 lire. Te 

I problemi economici riguardantii costi di trasmissione e di esercizio, quelli 
connessi all’organizzazione di un efficiente esercizio di ricezione e risposta, come 
pure uno studio più approfondito di tutto il problema sta nei programmi fu- 
turi del Centro. 


SUMMARY 


A description is given of the connections between memories, input-output units 
and the central part of the computer under construction at the Center of Pisa. The 
logical design takes into account the possibility of future enlargement of the external 
units and auxiliary memories of the Computer. The possibility has further been con- 
sidered of direct connections between the Center of Pisa and other scientific centers 
interested in the use of the computer. 
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IL NUOVO CIMENTO Vor. XII, N. 2 16 Aprile 1959 


Criteri elettronici 
di progettazione della calcolatrice del C.S.C.E. - I (*). 


G. CECCHINI (*), G. B. GERACE (**) e S. SIBANI (+%) 


Centro Studi Calcolatrici Elettroniche dell’ Università - Pisa 


(ricevuto il 28 Dicembre 1958) 


Riassunto. — Si espongono i criteri elettronici e costruttivi adottati 
- nella costruzione della « Macchina Ridotta » attualmente funzionante presso 
il Centro Studi Calcolatrici Elettroniche della Università di Pisa. 


La costruzione di questo calcolatore ha dato luogo ad una serie di decisioni 
di carattere sia generale (tipo, struttura), sia particolare (suddivisioni, elementi, 
tensioni e così via). 

Le prime decisioni, quelle che determinarono in sostanza il tipo di calco- 
latore, furono legate alle funzioni di tale macchina. Pensata come mezzo per 
risolvere problemi di carattere scientifico, dando quindi ad essa tutte le carat- 
teristiche di una macchina di tipo universale (general purpose), essa neces- 
sitava di alta velocità di calcolo e di grande flessibilità di programmazione. 

L'alta velocità di calcolo ha portato ad una macchina di tipo parallelo, 
nella quale cioè tutte le cifre binarie componenti il numero o l'istruzione ven- 
gono simultaneamente trattate. 

Si è inoltre voluta una semplicità concettuale di progetto, semplicità che si 
traduce in una grande semplificazione costruttiva. 

Perciò si è reso asincrono il funzionamento della macchina, adottando il 
criterio che le informazioni vengono trasmesse sotto forma di livelli e non di 


(*) Presentato. al XLIV Congresso della Società Italiana di Fisica, Palermo, 6-11 
Novembre 1958. 

(**) Distaccato dalla S.p.A Olivetti presso il C.S.C.E. 

(***) Distaccato dalla Sezione di Pisa dell’I.N.F.N. presso il C.S.C.E. 

(+) Attualmente presso la S.p.A. Olivetti. 
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impulsi; ciò ha portato ad una estrema semplificazione del controllo ed ha 
reso possibile la progettazione di un semplice addizionatore simultaneo ad alta 
velocità. 

Il parallelismo del calcolatore ha come conseguenza il forte numero di parti 
che si ripetono identicamente. Questo ha dato quindi luogo allo studio di una 
opportuna naturale suddivisione costruttiva del complesso, in modo da per- 
mettere una facile sostituibilità delle parti in caso di guasto, e, particolare 
non trascurabile, alla possibilità di una lavorazione, di un montaggio e di un 
controllo in serie, modo di operare che permette sopratutto di avere maggiori 
garanzie di sicurezza nella fase costruttiva. 

A questa concezione del calcolatore in parti staccate, ha fatto seguito il 
problema del collegamento di queste parti fra loro. 

Si è sistematicamente seguito il criterio di disaccoppiare gli organi fra loro, 
progettandoli in modo che presentassero sempre basse impedenze di uscita ed 
alte impedenze d’entrata. Si può in sostanza praticamente dire che nessun or- 
gano trasferisce corrente in un altro. Ne consegue che essendo state standar- 
dizzate le impedenze di entrata e di uscita, non esiste nessun problema per 
quanto riguarda l'accoppiamento dei vari organi fra loro, esclusione fatta per 
le capacità parassite dei fili di connessione. Va notato però che nei circuiti in 
continua le capacità parassite possono solo rallentare ma non pregiudicare il 
funzionamento della macchina. Infatti esse possono aumentare il tempo neces- 
sario per una determinata variazione di tensione senza alterarne però il valore 
finale. 

Inoltre i collegamenti fra organo ed organo sono stati realizzati con fili sem- 
plici perchè sono quelli che nelle condizioni di lavoro consentono le minori capa- 
cità parassite. D'altra parte l’uso di cavi coassiali chiusi sulla propria impe- 
denza caratteristica, avrebbe condotto ad erogazione di correnti ‘proibitive. 

La trasmissione delle informazioni sotto forma di livelli adottate per le grandi 
semplificazioni circuitali, presenta due aspetti negativi. 

Il primo è la maggiore dissipazione di energia del calcolatore ed il secondo 
è il problema dello spostamento dei livelli, fenomeno del resto ben conosciuto 
in tutte le trasmissioni in continua. 

Quest’ultimo problema da luogo ad un accurato studio sulla lunghezza 
massima delle catene di trasmissione, ossia sul massimo numero ammissibile 
di organi posti in serie in relazione all’elemento che utilizza i livelli stessi alla 
fine della catena. ; 

La soluzione di questo problema necessita una serie di decisioni sul tipo di 
tubo usato nella catena, sulla scelta delle tensioni e sulla loro tolleranza, sul 
tipo di resistore con riferimento alla sua tolleranza ed alla sua classe. 

Come si vede queste decisioni divengono un problema generale di macchina 
e possono venire prese solo quando si tenga conto di tutti i problemi. 

Nella scelta dei tubi si è sempre, ove possibile, osservato il criterio di usare 
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tubi a lunga vita, particolarmente studiati per il funzionamento a scatto (on-off). 

Ciò ha portato ad usare su larga scala il tubo E88CC di eccezionali presta- 
zioni, che presenta il vantaggio di potere erogare forte corrente con bassa ten- 
sione anodica e quindi con basse dissipazioni. Si noti il vantaggio dell’uso su 
larga scala di un solo tipo di tubo in modo da facilitare la manutenzione. 

Per alcune parti della macchina, specificatamente nella memoria, non si 
sono trovati in commercio tubi professionali specificatamente studiati per la- 
vorare on-off. Ci si è rivolti pertanto a tubi commerciali (ELS84) che, dato il 
modo di operare, ci garantivano analoga sicurezza di funzionamento nel tempo. 

Le tensioni di segnale, unificate in tutta la macchina, si ottengono mediante 
due tensioni di riferimento, il +10 ed il — 14. 

Queste tensioni sono asimmetriche, dato che la tendenza delle catene è di 
spostare più sensibilmente il livello inferiore verso l’alto di quanto avviene per 
l’altro livello. 

Le due tensioni, tenuto conto dei relativi spostamenti ora citati, sono state 
giudicate sicure per pilotare i circuiti a scatto e gli invertitori. 

Conseguenza di questa scelta è stata quella delle tensioni di alimentazione 
delle reti a diodi in +50 V — 50 V che consentono il migliore compromesso tra 
assorbimento di corrente e dissipazione della rete stessa. 

Queste due tensioni, oltre che per alimentare le reti a diodi, sono pure usate 
nei circuiti con tubi e permettono di minimizzare gli spostamenti dei livelli. 

Insieme al +50 V e — 50 V si sono utilizzati pure il +150 V e — 150 V. 

Si ottiene con queste due ultime tensioni il risultato di avere anodiche ade- 
guate per i tubi amplificatori ed un negativo convenientemente elevato al quale 
riportare i partitori di segnale senza troppo sacrificio per l’amplificazione. 

Per la memoria si è inoltre sentita la necessità di una nuova tensione di 
alimentazione dei tubi finali di potenza e per l’amplificatore di lettura, ten- 
sione fissata in -+200 V. 

Tutte le tensioni di alimentazione hanno una tolleranza del +2%; che corri- 
sponde ad un ragionevole compromesso tra le esigenze di stabilità delle tensioni 
ed il costo. 

Per ciò che si riferisce ai resistori si sono adottati elementi ad alta stabilità 
e con ristretta tolleranza per facilitare il problema degli spostamenti dei livelli. 

Il resistore è del tipo metallico depositato sotto vuoto. i 


SUMMARY 


The criteria adopted in planning the « Reduced Machine » at present working at 
the « Centro Studi Calcolatrici Elettroniche » of Pisa University are discussed. 


1037 


IL NUOVO CIMENTO Vor XIIT<N.#210, 16 Aprile 1959 


Criteri elettronici 
di progettazione della calcolatrice del C.S.C.E. - II (*). 


G. CECCHINI (*) e G. B. GERACE (**) 


Centro Studi Calcolatrici Elettroniche dell Università - Pisa 


(ricevuto il 28 Dicembre 1958) 


Riassunto. — Si discutono i risultati forniti dall’esperienza sulla « Mac- 
china Ridotta» e si accenna ad alcuni nuovi criteri adottati nella 
progettazione della « Macchina Definitiva » in costruzione presso il Centro 
Studi Calcolatrici Elettroniche dell’Università di Pisa. 


L’esperienza effettuata sulla « macchina ridotta » è stata di grande impor- 
tanza ai fini della progettazione elettrica della « macchina definitiva ». 

Essa ha permesso di apportare modifiche ad alcuni circuiti, di mutare cri- 
teri di scelta di componenti, di precisare problemi non del tutto chiariti, di dare 
al progettista garanzia e sicurezza su determinate scelte di impostazioni. 

Questa esperienza ci permette di pensare ora alla « macchina definitiva » 
come ad una macchina di sicure prestazioni e, pertanto, di comodo esercizio. 

Il lavoïo che ha condotto alla definizione del progetto della « macchina de- 
finitiva » è consistito sia in un riesame critico complessivo, nel senso sopra spe- 
cificato, sia nello studio di problemi nuovi, quali quello dell’analisi statistica 
del comportamento dei tubi con riferimento alle loro correnti di griglia e dell’in- 
troduzione di un comodo sistema di controllo della macchina. 


(*) Presentato al XLIV Congresso della Società Italiana di Fisica, Palermo, 6-11 
Novembre 1958. 

(”) Distaccato dalla S.p.A. Olivetti presso il C.S.C.E. 

(+) Distaccato dalla Sezione di Pisa dell’I.N.F.N. presso il C.S.C.E. 
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Gli argomenti trattati possono pertanto essere così elencati: 


1) Esame del comportamento dei componenti. 


2) Precisazione del concetto di durata nei tubi a lunga vita. 


3) Esame e modifica di circuiti dimostrativi critici. 
4) Analisi delle cause generatrici di disturbi di origine elettromagnetica. 


5) Analisi statistica dell'andamento delle correnti di griglia nei tubi in 
relazione all’invecchiamento. 


6) Introduzione di un sistema capace di controllare le diverse parti della 
macchina e di denunciare, con un determinato preavviso, eventuali disfunzioni. 


Passiamo ora ad una descrizione dettagliata dei vari punti. 


1) Tutti i componenti usati nella « macchina ridotta » hanno dato ottimi 
risultati, eccettuato i resistori. 

Nella scelta di questi così fondamentali componenti ci si era rivolti verso un 
tipo ottenuto mediante deposizione metallica sotto vuoto e garantito di qualità 
superiore alla classe 0.5 delle norme tedesche. 

In realtà si è osservato che dopo sole 1500 ore di vita (ovviamente a carico 
e temperatura nominali) circa il 20% di tali resistori è andato al di là del valore 
fissato dalla tolleranza ammissibile (questa tolleranza è data dal valore di tol- 
leranza iniziale più il valore della classe). 

Lo spostamento è unidirezionale e nel senso di accrescimento del valore. 

La spiegazione del fenomeno potrebbe probabilmente attribuirsi ad una 
imperfetta deposizione dello strato metallico, causa di un sovrariscaldamento 
di alcune zone e conseguente evaporazione dello strato. 

Questo modo di vedere potrebbe venire convalidato dal fatto che le devia- 
zioni più alte si hanno per resistori di valore più grande. 

Il componente è stato pertanto sostituito con un tipo di precisione a pellicola 
di carbonio lungamente stagionata, che, da prove sperimentali effettuate presso 
il nostro laboratorio, ha mostrato un ottimo comportamento. 


2) Nella progettazione iniziale della « macchina ridotta » non si era se- 
guito un concetto unitario che tenesse conto delle variazioni delle caratteri- 
stiche dei tubi durante il funzionamento. 

Era stato dato ad ogni buon conto, un margine di sicurezza, ma senza un 
approfondito esame del comportamento dei tubi. Inoltre questo margine non era 
stato unificato, ma era stata lasciata la decisione al buonsenso del progettista. 
Già prima dell’inizio della costruzione della « macchina ridotta » si prese in 
esame questo problema e si esaminarono i circuiti nel loro comportamento con 
variazioni percentuali nelle correnti, da inizio a fine vita, da +15% a — 30%; 
limiti che risultarono essere stati, nella quasi totalità dei casi, notevolmente 
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superati. Dopo la costruzione della « macchina ridotta » il problema fu decisa- 
mente affrontato e fu messo in luce il fatto che un tubo viene considerato a 
fine vita quando ha raggiunto il — 30% della G, nominale, corrispondente, 
ad esempio per il tipo E88 CC che rappresenta la quasi totalità dei tubi usati 
e per le condizioni di lavoro normalmente usate, ad un — 50% del valore no- 
minale delle correnti. | 

Per la macchina definitiva questo è stato considerato il limite inferiore del 
tubo E88 CC, per cui cioè i diversi circuiti possono ancora lavorare con sicurezza. 
Analoghe considerazioni si estendono ‘per tutti gli altri tipi che però, ripetiamo, 
corrispondono ad una ‘piccola percentuale del totale dei tubi. 


3) L'osservazione del comportamento della « macchina ridotta » ha por- 
tato alla determinazione di circuiti dimostratisi critici o non adeguatamente 
progettati. Si sono pertanto modificati nel nuovo progetto il registro standard, 
il distributore di impulsi e quello dei segnali, e l’invertitore standard. 

Per il resto si è rifatta la verifica con i nuovi limiti precisati sopra, tenendo 
opportunamente conto del fenomeno del quale si parlerà al punto 5) e riscon- 
trando l’ottimo proporzionamento dei circuiti. 


4) AlVingresso dei decodificatori della Memoria si sono riscontrati disturbi 
di entità tale da indurre in errore il calcolatore. Questi decodificatori sono col- 
legati con l’unità di calcolo mediante fili di lunghezza di circa 2 metri posti in 
vicinanza di altri fili portatori di informazioni. Si sono effettuate lunghe prove 
che hanno definito di origine elettromagnetica la causa del fenomeno. 

Le brusche variazioni di livello nei fili portatori di informazioni producono 
variazioni di corrente che, trascurando la componente induttiva, sono funzione, 
all’istante di inizio dal fenomeno, solo della pendenza del fronte di tensione dato 
dal generatore e della sua resistenza interna. 


5) Si è sperimentalmente dimostrato che le correnti di griglia del tubo 
E88CC hanno una variazione con l’invecchiamento del tubo, nel senso di diminuire 
allorquando diminuiscono le correnti di placca del tubo stesso. Questo problema 
quindi va ricollegato con il fenomeno della variazione di emissione del catodo 
col tempo, come lo è per le correnti anodiche. 

Si sono fatte lunghe ed accurate prove per accertare le entità di queste va- 
riazioni e si è notato lo stretto collegamento esistente fra le due variazioni di 
corrente di griglia e placca. Questo fenomeno è stato utilizzato per diminuire 
l’effetto delle variazioni delle caratteristiche dei tubi, mediante l’effetto di 
compensazione esercitato da una resistenza posta in serie alla griglia stessa. 

È stato reso pertanto più agevole il problema della minimizzazione degli spo- 
stamenti dei livelli positivi di segnale lungo la catena di trasmissione in continua 
ed in particolare si è resa possibile una più larga stabilizzazione delle correnti 
negli amplificatori funzionanti con griglia positiva (invertitori ed elementi à 
scatto). 
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6) Si è infine studiato un semplice, efficace sistema di controllo di macchina 
che permetterà di seguire l'andamento complessivo nel tempo di ogni catena 
di trasmissione e di ogni organo e di individuare l’elemento che si avvicina alla 
zona di criticità di operazione con largo margine di tempo, in modo da facili- 
tare al massimo l’operazione di sostituzione. 

Questo sistema verrà corredato da programmi appositamente studiati che 
permettono di realizzare le condizioni più sfavorevoli per l’organo da control- 
lare. Il sistema si basa sul concetto di disturbare di volta in volta le varie 
catene nelle quali la macchina può essere suddivisa e di determinare l’entità 
del disturbo necessario per indurre in errore il calcolatore. 

La determinazione del punto di errore può essere fatta con sicurezza solo 
mediante un opportuno programma che permetta alla macchina stessa di av- 
vertire l’operatore quale è l’elemento singolo della catena in esame che ha pro- 
vocato l’errore. Prove di controllo fatte sulla « macchina ridotta » hanno dato 
ottimi risultati e lasciano sperare in una semplice, comoda e veloce routine di 
controllo giornaliero della « macchina definitiva ». 


SUMMARY 


Experience acquired from the « Reduced Machine » is described together with some 
new criteria adopted in planning the « Final Machine » which is being built at the 
« Centro Studi Calcolatrici Elettroniche » of Pisa University. 


9 - Il Nuovo Cimento. 
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Dispositivo atto a selezionare nella radiazione cosmica 
eventuali particelle di massa intorno a 550 m, (*). 


M. Conversi, G. M. De’ MUNARI, A. EGIDI, E. FIORINI, S. RATTI, C. RUBBIA 


C. Succt e G. TORELLI 


Istituti di Fisica delle Università - Milano e Pisa 
Istituto Nazionale di Fisica Nucleare - Sezioni di Milano e Pisa 


(ricevuto il 28 Dicembre 1958) 


Riassunto. — Si descrive il dispositivo utilizzato al fine di controllare 
l’esistenza recentemente segnalata da ALTHANIAN et al. (1) di particelle di 
massa intorno a 550 m,. Si discutono le caratteristiche dell'apparato le 
cui prestazioni sono tali da garantire una completa eliminazione del 
«fondo » ed un sicuro riconoscimento delle eventuali particelle ricercate. 


1. — Una comunicazione di ALIHANIAN et al. al Congresso di Mosca del 
1956 (1) ha risvegliato l’interesse sul problema delle particelle di massa inter- 
media tra la massa del pione e quella del mesone K, l’esistenza delle quali si 
sospettava in base ad indicazioni non conclusive di vari autori (2). ALTHANIAN 
edi suoi collaboratori, riportando i risultati di una serie di esperienze sullo spettro 
di massa delle particelle cariche della radiazione cosmica, hanno presentato 


(*) Presentato al XLIV Congresso della Società Italiana di Fisica, Palermo, 6-11 
Novembre 1958. 

(1) A. I. ALIHANIAN et al.: Intern. Conference on High Energy Physics (Mosca, 
Maggio 1956). Cfr. A. I. ALIHANIAN, N. V. SosTAKOvIŸ, A. T. DADAIAN, N. V. FE- 
DOROV e B. N. DERIAGIN: Zurn. Éksp. Teor. Fiz. (U.S.S.R.), 31, 955 (1956); Soviet 
Phys. Journ. Exp. Theor. Phys., 4, 817 (1957) e bibliografia ivi contenuta. 

(2) J. G. RETALLACK and R. B. BropE: Phys. Rev., 75, 1716 (1949); R. B. BRODE: 
Rev. Mod. Phys., 21, 37 (1949); P. R. DanieL and H. D. PERKINS: Proc. Roy. Soc., 
A 221, 351 (1954) e bibliografia ivi riportata. 
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in effetti una concreta evidenza in favore della esistenza di particelle di massa 
prossima a 550 masse elettroniche (*). 

11 dispositivo sperimentale impiegato da ALIHANIAN et al., è costituito da 
uno spettrometro magnetico interposto tra due camere di Wilson a lastre. La 
misura della massa viene ricavata dalla curvatura della traiettoria della par- 
ticella e dal «range » misurato nella camera inferiore che contiene lastre di 
piombo. La camera a lastre superiore consente di riconoscere se le particelle 
in esame provengono o meno da interazioni locali. Lo spettro di masse ottenuto 
per le particelle che non provengono da interazioni locali presenta un picco 
molto netto, costituito da 11 particelle con masse distribuite intorno a 530 mo 

ALIHANIAN et al. concludono in definitiva che esistono nella radiazione 
cosmica particelle dotate delle seguenti caratteristiche: 


a) massa prossima a 550 m,; 
b) carica prevalentemente negativa (9 particelle negative e 2 positive); 


c) produzione non locale; (e quindi vita media presumibilmente abba- 
stanza elevata); 


d) percentuale relativa ai mesoni u nello stesso intervallo di «range » 
~ 0.5%. 


La bassa frequenza delle particelle 550 e la mancanza di informazioni sul 
decadimento, ci hanno indotto a dare particolare peso, nella progettazione del 
nostro apparato, al problema della eliminazione del « fondo ». Quest’ultimo è 
essenzialmente dovuto a possibili interazioni nucleari delle particelle selezionate 
dall’anparato. Più precisamente se, come nell’esperienza del gruppo di Alihanian, 
le determinazioni di massa sono ricavate da misure di percorso e di curvatura 
in campo magnetico, per le particelle che subiscono interazioni nucleari si trova 
una massa apparente più grande di quella reale (3). Con questo sistema di misura 
il pericolo di registrare particelle 550 spurie è rappresentato quindi da mesoni x. 
Se la determinazione di massa si basa invece su misure di ionizzazione e di 
percorso residuo, le particelle che subiscono interazioni presentano una massa 
apparente inferiore a quella reale (*). In tal caso il pericolo di registrare parti- 
celle 550 spurie deriva essenzialmente da protoni interagenti. 

Nella esperienza di cui in questa comunicazione si descrive il dispositivo, 
abbiamo scelto questo secondo metodo di misura per ottenere un angolo solido 


(*) Indicate nel seguito con l’abbreviazione « particelle 550 ». 

(3) M. Conversi, G. M. De’ MUNARI, A. Earp, E. Fiorini, S. RATTI, C. Russia, 
C. Succi e G. TORELLI: Ricerca sull’esistenza di particelle di massa intorno a 550 m,, 
nella radiazione cosmica. Comunicazione al XLIV Congresso della S.I.F. Cfr. anche 
C. RuBBIA: Tesi di laurea (Università di Pisa, 1957). 
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di accettazione sufficientemente grande. Tale esperienza (che si sta ora conclu- 
dendo al Laboratorio Raggi Cosmici del Sabbione (Novara) a 2550 m s.m.) è 
stata intrapresa al fine di controllare l’esistenza delle particelle 550 ed eventual- 
mente di studiarne le proprietà. 


2.:- Il dispositivo sperimentale impiegato in una prima serie .di misure è 
presentato in Fig. 1: A, B, C, D ed E sono contatori di Geiger-Müller; 8, ed 


Supporto 
di legno 


id Scintillatore Sc 


ni 


TT] 
Ti 
| 


Di 


À 


Scintillatore ri anticoincidenza Sa 


Fig. 1. — Disposizione dei rivelatori e delle camere di Wilson nella 1° serie di misure 
(Q = 10.5 em? sr). 


$, scintillatori liquidi; ©, un contatore di Cerenkov; CI e CR due camere di 
Wilson di cui la seconda a lastre. Le coincidenze 0, D,, 0, Di, 03 D, tra i gruppi 
di contatori C e D indicati in figura, definiscono l’angolo solido totale di accet- 
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‘tazione, interamente coperto, al di sotto della camera di Wilson inferiore (CR), 
dallo scintillatore $, posto in anticoincidenza. In tal modo si selezionano le 
particelle che terminano il loro percorso nella camera CR dopo aver attraver- 
sato la camera superiore CI. 

Il contatore di Cerenkov C,, in anticoincidenza, impedisce che la compo- 
nente molle ed i mesoni leggeri (7 e u) comandino le camere. Alla eliminazione 
della componente molle contribuiscono anche i contatori Æ in anticoincidenza 
ed i contatori A e B in coincidenza, separati da uno schermo di Pb di 1.5 cm 
di spessore. 

La determinazione di massa delle particelle che si arrestano nelle lastre di 
CR viene generalmente (*) ricavata dalla misura della perdita di energia nello 
scintillatore proporzionale S. e dalla misura del percorso residuo in CR. 

Data la buona risposta dello scintillatore S, la camera CI, inizialmente 
destinata a misure di ionizzazione, è stata usata solo per controllare l’allinea- 
mento delle tracce delle particelle in esame nelle due camere onde facilitare 
il riconoscimento delle eventuali interazioni nucleari tra lo scintillatore $, e la 
camera CR. 

Un ulteriore forte criterio per la identificazione delle particelle in esame è 
fornito dalla camera CR che contiene 17 lastre di plexiglas di 0.3 cm di spes- 
sore. In effetti: 


a) il sottile spessore delle lastre facilita il riconoscimento delle interazioni; 


b) l’elevato numero dei «gaps» a disposizione permette una stima ab- 
bastanza accurata della massa in base alla variazione di ionizzazione ed allo 
scattering della particella (3); 


e) il basso numero atomico del materiale usato per le lastre rende minimo 
lo scattering coulombiano e quindi più evidenti le deflessioni dovute alle inte- 
razioni nucleari. 


Per una seconda serie di misure attualmente in via di conclusione, il dispo- 
sitivo sperimentale precedentemente impiegato è stato modificato come in Fig. 2 
allo scopo di ottenere un maggiore angolo solido di accettazione (17.5 cm? ster 
in luogo di 10.5 cm? sr), ed una ulteriore informazione sulla massa delle parti- 
celle osservate. Come risulta dalla figura, la camera di Wilson CI è stata sosti- 
tuita da un secondo contatore di Cerenkov €, il cui impulso non viene utiliz- 
zato per la selezione elettronica del comando bensì in sede di analisi. Il 
contatore 0, permette un immediato riconoscimento degli eventi generati da 
protoni. Infatti lo spessore che lo separa dalla camera CR è calcolato in modo 
che le particelle che si arrestano in CR perdendo energia solo per ionizza- 
zione (cioè senza subire interazioni nucleari) eccitano €. soltanto se la loro 
massa è inferiore a — 1000 m,. D'altra parte un protone che produce luce 
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attraversando C,, perde nel contatore proporzionale S, una energia decisa 
mente inferiore a quella di una particella 550 che si fermi nella camera CR. 
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Fig. 2. — Disposizione dei rivelatori e della camera di Wilson nella 22 serie di misure 
(OF s7-o) Cm 182). 


Gli spessori compresi in entrambe le serie di misure tra ©, e CR sono 
tali da escludere che una particella 550 che si fermi in CR produca luce in C,. 
È opportuno osservare, a questo proposito, che è stato necessario aggiungere, 
nella seconda serie di misure, 3 cm di Al al di sotto del contatore ©, per otte- 
nere una completa eliminazione dei mesoni y. 

Nelle misure della seconda serie, quindi, i due contatori di Cerenkov ed i 
due scintillatori costituiscono un selezionatore di masse che accetta soltanto 
particelle di massa intermedia tra i mesoni leggeri e i protoni. 
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Le caratteristiche dei rivelatori S., S,, ©: e C, sono riportati in Tab. I; 
in Tab. II sono riportate le principali caratteristiche delle camere. 


TABELLA I. — Caratteristiche dei rivelatori a scintillazione ed a luce di Cerenkov usati 
nell’apparato di Fig. 2. 
Rive- |Volume sen- x Materiale | Wave | Diffu- : | Rendi- | 
latore | sibile (cm8) Tipo sensibile | Shifter sori Fototubi mento 
2 6 x 6292 0.998 
(DE 62x 24 12| Cerenkov | Acqua Nessuno| MgO Du Mont | (g> 0.8) | 
Ss 12 x 6291 | 0.998 
+ 
C, (*) | 62X24X 18 | Cerenkov | Acqua Nessuno; MgO Du Mont | (8> 0.8) | 
scintil Toluolo oe 
8S, | 48x24x5 ioe a +3.5g/ |xNPO | Al 4 x 629 Dai 
: atore Pacatert..| Du Mont | 
Scintil Toluolo 4 
Sa |106x62x6| °° | 35 el la NPO | Al 6x6292 | 0.997 
latore | Du Mont | 
Paraterf. 
(*) Usato come contatore di anticoincidenza (C,) nella disposizione di Fig. 1. 


TABELLA II. — Caratteristiche delle camere di Wilson. 


Camera di Wilson 


CI 


CR 


Dimensioni utili em* 
Parete superiore plexiglas 
Parete inferiore 
Lastre interne 


50 x 40 x 12.5 

0.8 cm 

1.8 cm plexiglas 

4 verticali da mm 3 
di Al 


110 x 110 x 50 

0.8 em 

2.5 cm Al 

17 orizzontali da 
mm 3 di plexiglas 


8. — Lo schema di principio dell’apparato elettronico è mostrato in Fig. 3. 
Gli impulsi dei contatori C., 0, ed S, sono inviati ad altrettanti « gates » aperti, 
per la durata: di circa 100 ns, dagli impulsi del contatore S,. In tal modo è 
possibile evitare l’affollamento di impulsi che diversamente si avrebbe sui cir- 
cuiti di anticoincidenza e sull’amplificatore A. 

L'energia persa dalle particelle attraversando lo scintillatore S, viene regi- 
strata fotografando l’impulso di questo rivelatore (opportunamente amplificato) 
sul tubo di un oscilloscopio; la «sweep » dell’oscilloscopio è comandata dallo 
stesso impulso che produce l’espansione delle camere di Wilson (Fig. 3). Sulla 
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Cerenkov C 


Impulsi di 


coincid. (AB) ” 


Cerenkov Ca 777] 


Asse y 
Tektronix 515 
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ABCDS, (Ca*Cq) 


| 


Comando camere 


Scintillatore S> 


Impulsi di 
coincid. (C D)= 


(C,D,*C,D,*CD,) 


Scintillatore S, 


(a) — Canale di anticoincidenza 
A — Amplificatore lineare 
Ant ~ Anticoincidenza 

C = Coincidenza 

F = Formatore di impulsi 

G > Gate, 

M + Mescolatore di impulsi 
R — Ritardo 


Fig. 3. - Schema a blocchi del dispositivo elettronico che seleziona gli impulsi di comando. 


stessa «sweep », nella seconda serie di misure, viene inviato anche l’impulso 
ritardato del contatore C,. 
Il dispositivo elettronico è stato progettato ponendo particolare attenzione 


120F | 


110} 
100} 


3 
go! > 
Ss 


nm 1evento 


a Hlmin 


=> 
aL 
rok 


al problema della stabilità, essenziale in un 
apparato destinato a lunghi periodi di fun- 
zionamento. Particolarmente critici, perchè 
da essi dipende l’esattezza della misura di 
massa, sono il rendimento luminoso dello 
scintillatore S, e la stabilità dell’amplifi- 
catore ad esso connesso: per queste ragioni 
è stato giornalmente eontrollato lo spettro 
di ampiezza degli impulsi di S, prodotti 
dalla radiazione cosmica incidente nell’an- 
golo solido di accettazione. Tali spettri 
(Fig. 4) permettono anche di calcolare 
volta per volta l’errore della valutazione 
della perdita di energia nel contatore S,. 

Tutti gli altri controlli, utili per accer- 
tare le perfette condizioni di lavoro delVap- 
parato, si possono eseguire senza pertur- 
barne il funzionamento. In particolare il 
controllo dei contatori GM viene eseguito 
a distanza con il circuito schematicamente 


Fig. 4. - Esempio di spettro di controllo dello 
scintillatore S,. 
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indicato nella Fig. 5. Il controllo del rapporto di espansione della camera di 
Wilson è effettuato con una macchina fotografica tale da permettere il rapido 


esame del fotogramma e quindi la correzione del rapporto prima dell’espan - 
sione successiva. 


+A.T. +150 
Canale di 
22M _ Mescolazione 
} 
Uscita di 
mescolazione 
| 
| = 
+150 
ee] 
K ; a 
+50 Canale di Uscita di 
controllo controllo 
1M 
e 


lb 


Fig. 5. — Schema del circuito per il controllo a distanza dei contatori di Geiger-Miiller. 


È opportuno sottolineare, in conclusione, che il dispositivo descritto, pur 
non consentendo accurate determinazioni di massa, è in grado di fornirci una 
quantità di informazioni tale da garantire una completa eliminazione del fondo 
ed un riconoscimento sicuro delle eventuali particelle 550. 


SUMMARY 


The apparatus employed in a search for the 550-mass particles reported by ALIHA- 
NIAN et al. (1) is described. Unwanted particles are eliminated electronically by a 
suitable arrangement of scintillation — Cerenkov —, and GM-counters. The recorded 
particles are observed in a large multiplate Wilson chamber. The mass is essentially 
deduced by the combined measurements of range in this chamber and dE/dx in a 
proportional liquid scintillator. It is shown that the background of spurious 550-mass 
particles is entirely negligible and that genuine 550-mass particles could not be missed. 
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Impiego della « macchina ridotta» del C.S.C.E. di Pisa 
nella soluzione di alcuni problemi (*). 


E. FABRI (**) e L. GUERRI (**) 


Centro Studi Calcolatrici Elettroniche dell Universita - Pisa 


(ricevuto il 28 Dicembre 1958) 


Riassunto. — Si descrivono i calcoli effettuati o in corso di elaborazione 
per la «macchina ridotta » e si riferisce su alcune esperienze di moderni 
metodi di calcolo numerico. 


1. — Introduzione. 


La Macchina ridotta del C.S.C.E. è entrata in funzione in modo completo 
verso la metà del 1957. Da allora fino alla fine del 1958 essa è stata disponibile 
per calcoli per una larga parte del tempo. Sebbene la macchina non fosse stata 
originariamente prevista per un impiego effettivo, ma solo per collaudi e prove 
tecniche, è stato possibile svolgervi un discreto lavoro, anche se le limitazioni 
cui si farà cenno in seguito hanno pesato sulla possibilità di eseguire calcoli 
di una certa complessità. Esporremo ora, seguendo all’incirca un ordine cro- 
nologico, il lavoro svolto sulla Macchina Ridotta (M.R.) indicando per ogni 
calcolo origine, difficoltà, tempo richiesto per la programmazione e per l’ese- 
cuzione, e ogni altro dato che possa riuscire interessante. 


(*) Presentato al XLIV Congresso della Società Italiana di Fisica, Palermo, 
6-11 Novembre 1958. 

(**) Attualmente all’Istituto di Fisica dell’Università di Pisa. 

(*»*) Distaccato dalla Sezione di Pisa dell’I.N.F.N. presso il C.8.C.E. 
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2. Studio di strutture di cristalli. 


Questo studio si riconduce ad un’analisi diretta e inversa di Fourier. 
Nel caso dell’analisi diretta si doveva calcolare: 


N N 


F(a, y) = 2 DIL rs COS 2r(rx + sy), 


assegnati gli a,, e N (con N<15) e: 


a=0d, 1d, 2d.. 50d, 
y=04, 1d, 2d.. 100d, 


4 — 0.02 oppure d = 0.01. 

I valori distinti del coseno sono 100 (per d = 0.02 ne servono solo 50). Questi 
100 valori del coseno sono stati incorporati nel programma. 

Per d = 0.02 si ha un reticolo di 1326 punti. L'intero calcolo comporta 
(N +1)(2N +1)-1326 moltiplicazioni. La moltiplicazione, non essendo disponi- 
bile come istruzione, era eseguita come sottoprogramma. Tempo di calcolo: 
1 ora circa di cui 4 impegnato per la stampa dei risultati fatta, carattere per 
carattere, a mezzo di telescrivente. 

Nel caso d= 0.01 le cifre e i tempi sopraindicati vanno moltiplicati per 4. 

Per l’analisi inversa si doveva calcolare: 


N 
G(h, k) = > a,(h, k) cos 2x(hr, + ky.) , 
r=0 
assegnati: 
a) N(con N<20); 
b) l'insieme dei valori di h e k e dei relativi coefficienti a,(h, k) (r = 0, 
MZ EN) = 
cha © Y, (M=0, 1, 2... N) compresi fra 0 e 1. 
Il coseno era calcolato per mezzo di un sottoprogramma. 
Tempo di calcolo: variabile fra 5 e 10 min. 


3. — Analisi periodale di frequenze critiche ionosferiche e correlazioni con i 
dati sulle attività solari. 


Questo è il calcolo che ha-occupato per più lungo tempo la macchina causa 
la notevole massa di dati da elaborare. Il calcolo constava di due parti: la 
prima era l’analisi periodale delle frequenze critiche ionosferiche, con la ricerca 
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delle componenti a periodo 12h e 24h nonchè delle componenti a periodo annuale 
e dell'andamento secolare. Da un punto di vista matematico il problema non 
presentava alcuna difficoltà; c’era solo la necessità di introdurre un gran nu- 
mero di dati (circa 25000 numeri di tre cifre) e di stampare un numero ancora 
più grande (circa 80000 di quattro cifre). Tutto ik-ealcolo è stato ripetuto due 
volte per controllo, e tre o quattro volte nei punti dove si scoprivano errori. 

Il tempo effettivo di calcolo era trascurabile rispetto a quello di stampa, 
nonostante che per ogni numero stampato occorressero in media quattro mol- 
tiplicazioni più un numero imprecisato di addizioni e sottrazioni. Questa prima 
parte del calcolo ha richiesto in totale circa 60 ore di macchina, distribuite 
però in parecchi giorni a causa di errori, interruzioni varie, necessità di con- 
trolli, ecc. Il lavoro di programmazione può essere valutato a 15 giorni di due 
persone. 

La seconda parte del calcolo era matematicamente più semplice della pre- 
cedente, ma presentava difficoltà di organizzazione. Occorreva infatti utiliz- 
zare i risultati della prima parte (già disponibili su nastro perforato) per il 
:alcolo delle somme e somme di quadrati per ogni singola serie di dati; occor- 
reva poi calcolare le somme dei prodotti dei termini di ciascuna delle serie con 
quelli di ciascuna serie di dati dell’attività solare (preventivamente preparati 
su nastro). Non è il caso di entrare ora in dettagli sulla questione: basta accen- 
nare che il nastro passato attraverso l’organo di entrata della macchina si mi- 
sura in km. I risultati da stampare erano circa 15000 numeri di dieci cifre e il 
tempo impiegato è circa lo stesso della prima parte, con le stesse avvertenze. 
Il lavoro di programmazione è stato invece alquanto minore: circa 10 giorni 
di una persona. i 


4. — Calcolo delle autofunzioni e dei livelli fondamentali di sistemi atomici. 


Il lavoro in questione è stato studiato e programmato a titolo di prova da 
un gruppo dell’Istituto di Chimica Fisica, nel caso della molecola-ione di idro- 
geno. Il calcolo era condotto con un metodo variazionale-iterativo e ha dato. 
risultati soddisfacenti. Non siamo in grado di valutare il tempo di programma- 
zione, che non si discosta però come ordine di grandezza da quelli già indicati; 
il tempo di calcolo era molto breve, di alcuni secondi nel caso effettivamente 
eseguito. 


5. — Calcolo dell’equilibrio di un reticolo cristallino. 


Il problema era riducibile alla ricerca del minimo di una funzione algebrica 
di 4 variabili, alquanto complicata causa la presenza di molte radici quadrate. 
In questo caso non esistevano questioni derivanti dalla mole del lavoro, ma le 
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limitazioni principali sono derivate dalla scarsa precisione ottenibile (l’uso 
delle tecniche di « doppia precisione » era impraticabile perchè avrebbe com- 
plicato il programma in modo eccessivo rispetto alla importanza del lavoro). 
Il programma ha effettivamente funzionato, ma senza fornire risultati utili 
giacchè le variazioni dovute ad arrotondamenti ed approssimazioni nasconde- 
ivano il minimo reale della funzione. Il tempo di macchina era in questo caso 
trascurabile; la programmazione ha richiesto 10 giorni di due persone. 


6. — Aggiustamento su dati sperimentali di una funzione di prova. 


Il calcolo, condotto con il metodo dei minimi quadrati, utilizzava per la 
ricerca del minimo lo stesso programma impiegato nel caso precedente (cal- 
‘colo dell’equilibrio di un reticolo cristallino). Il programma ha fornito i risul- 
tati voluti, che hanno però mostrato che la scelta della funzione di prova era 
inadeguata. Il tempo di macchina è stato di qualche ora, quello di programma- 
zione di circa una settimana per una ‘persona. 


7. — Ricerche sui metodi di Monte Carlo. 


Sono state condotte varie ricerche sull’impiego di tecniche probabilistiche 
nella risoluzione dei vari problemi. Il lavoro è stato orientato in primo luogo 
verso la ricerca del più conveniente metodo di produzione di numeri « pseudo- 
casuali » (*). Si è trovato che buoni risultati si ottenevano con il metodo noto 
.col nome di « serie di Fibonacci », purchè si lavorasse in doppia precisione. 
Ciò per evitare il presentarsi di quasi-periodi troppo brevi nella successione, 
che avrebbero impedito la costruzione di statistiche abbastanza ampie senza 
errori sistematici. 

Col metodo accennato, e con altri sperimentati, è stata programmata la 
risoluzione di sistemi di equazioni lineari. 


(*) Conviene accennare che l’impiego di un metodo di Monte Carlo richiede la pro- 
duzione di eventi casuali, che si può sempre ridurre a quella di numeri casuali. Senza 
ricorrere a effettive estrazioni o sorteggi, si può ricorrere a tavole di numeri casuali 
già esistenti; tale procedimento non è però indicato per una macchina con memoria 
piccola ed entrata non molto veloce, come la nostra M.R. L’altra possibilità è quella 
di far produrre alla macchina, secondo un processo deterministico, successioni di nu- 
meri che abbiano « apparenza casuale », nel senso che la regolarità della successione 
sia completamente nascosta sotto un’apparente irregolarità. I numeri così ottenuti 
si dicono pseudo-casuali. Perchè lo si possa accettare, un metodo di produzione di 
numeri pseudo-casuali deve superare appositi controlli, sui quali non ci possiamo qui 
«soffermare. i È 
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Il risultato è stato il seguente: a causa delle limitazioni della macchina non 
era ‘possibile rendere il metodo di Monte Carlo competitivo con quelli alge- 
brici; ciò a causa dell’eccessivo tempo necessario a produrre il singolo evento 
casuale. 

Un dato indicativo è il seguente: la risoluzione all’1% di un sistema di 10: 
equazioni si otteneva in alcuni minuti, ma-il tempo necessario cresceva col 
quadrato della precisione richiesta. Il metodo di Monte Carlo è conveniente 
per precisioni non troppo elevate e con molte equazioni (il tempo di calcolo 
va infatti col quadrato del numero di equazioni, mentre coi metodi algebrici 
va col cubo), ma la M.R. non consentiva, causa la scarsità di memoria, di spe- 
rimentare i casi interessanti. 

Strettamerte legata alla soluzione dei sistemi lineari è l’inversione di matrici, 
che è stata tentata con risultati del tutto analoghi. 

Altre prove sono state condotte su equazioni integrali (problemi di soprav- 
vivenza) e su problemi direttamente probabilistici, sempre con risultati simili: 
si ottenevano facilmente delle soluzioni grossolane, ma il tempo di calcolo di- 
veniva proibitivo se si cercava di raffinare i risultati. 

Un'altra difficoltà pratica che conviene accennare è risultata dall’even- 
tualità di errori della macchina: mentre in un calcolo ripetitivo un errore può 
influire su uno solo o pochi dei molti risultati finali, nel caso dei metodi di Monte 
Carlo è necessario che nessun errore si verifichi per tutto il tempo del calcolo, 
pena la perdita totale del lavoro già fatto. 


8. — Determinazione esatta delle autofunzioni del momento angolare di più 
elettroni. 


Questo calcolo è descritto in dettaglio in altra comunicazione (*), alla quale. 
si rimanda per ogni notizia. 


9. — Varie. 


È appena il caso di accennare al lavoro di programmazione preliminare 
consistente nella elaborazione di tutti i sottoprogrammi che sono stati impiegati 
nei calcoli descritti; va solo detto che il relativo tempo di programmazione 
— abbastanza considerevole, data anche l’inesperienza iniziale — non è stato 
considerato nelle stime date sopra. 


(*) E. ABATE ed E. FABRI: Impiego di una calcolatrice elettronica nella ricerca delle: 
autofunzioni del momento angolare orbit ale in accoppiamento di Russel-Saunders. 
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Altro tempo è stato preso dalla compilazione di programmi di prova e dia- 
gnostici, sia pure molto semplici, necessari per controllare la macchina e sco- 
prire eventuali guasti. 

Altri programmi effettivamente eseguiti sulla M.R. vanno meglio catalo- 
gati come esercizi e studi di programmazione che come applicazioni pratiche; 
ad esempio lavori di teoria dei numeri (divisori di un numero, M.C.D., numeri 
‘primi, numeri perfetti ecc.) o questioni a carattere combinatorio (permuta- 
zioni ecc.). 

Tra questi ultimi programmi merita di essere segnalato per la sua comples- 
sità logica e per un certo interesse teorico quello per la costruzione di tautologie 
del calcolo delle proposizioni. Si tratta di un problema di logica simbolica, im- 
mediatamente riducibile ad un problema combinatorio: disporre dei simboli 
di proposizioni, dei simboli di connettivi logici e delle parentesi in modo che 
le espressioni risultanti, se interpretate secondo opportune regole, risultino for- 
mule di proposizioni identicamente vere (tautologie). Il programma ha potuto 
essere fatto in pochi giorni, ma il tempo di macchina necessario per ottenere 
tautologie non banali è risultato eccessivamente lungo (giorni o mesi). 


SUMMARY 


A report is given on problems solved or being solved on the «reduced machine » 
and on some experiments with modern computing techniques. 
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Un nuovo tipo di discriminatore differenziale di impulsi rapidi (’). 


C. RUBBIA e G. TORELLI 


Istituto di Fisica dell’ Università - Pisa 
Istituto Nazionale di Fisica Nucleare - Sezione di Pisa 


(ricevuto il 28 Dicembre 1958) 


Riassunto. — Si descrive un nuovo tipo di discriminatore differenziale 
capace di funzionare con elevate frequenze di ripetizione nonostante la 
notevole semplicità della sua struttura. 


È in corso di realizzazione un discriminatore differenziale multicanale di 
nuovo tipo che rispetto agli strumenti del genere oggi esistenti, presenta note- 
voli vantaggi per semplicità ed economia, pur potendo 
giungere a frequenze di registrazione dello stesso ordine 
di grandezza. Esso è basato sull'impiego di un tubo 
Philips a deviazione di fascio (il tubo « E 80 T » costruito 
come confrontatore di frequenza per apparati televisivi). 
Questo ha la struttura di un tubo a raggi catodici 
(Fig. 1) con fascio laminare, mancante naturalmente di 
una delle coppie di placchette deflettrici. Lo schermo è 
costituito da una placca che riceve il fascio di elettroni 
solo quando la differenza di potenziale tra le due plac- 
kg? f fF chette D e D' è inferiore a AV, (Fig. 2). Forzando la fo- 

calizzazione e diminuendo la tensione di schermo si 
può ridurre la larghezza della curva di risposta senza 
preoccupanti riduzioni della corrente massima del 
tubo (Fig. 3). 

L’utilizzazione di questo tubo come discriminatore differenziale è relativa- 
mente semplice; si blocca una delle placchette, per esempio la D', al normale 


a €80T 


Fig. 1. — Struttura del 
tubo E 80 T. 


(*) Presentato al Congresso della Società Italiana di Fisica, Palermo, 6-11 No- 
vembre 1958. 
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potenziale di funzionamento V,, e si richiama l’altra placchetta, la D, ad un 
potenziale V, = V,- V ove la tensione V è arbitraria. Se sulla placchetta D 
si invia un impulso rettangolare di ampiezza V proporzionale all'impulso in 
esame, la E 80 T viene a trovarsi in condizioni di conduzione solo se l'ampiezza 


x 


dell’impulso è compresa tra V—AV e V+ AV. 


Y, =100V 


Vg = 70V 
Vig? gh) =250V 


Z| \ 
20 40 60 80 100 


120 


Fig. 2. — Caratteristiche del tubo E 80 T 


fornite dalla Philips. 


(mA) |Va =170V 
12/4 vy? = 80V 
1 Vigg4y=170V 

=22V 
pale = 


On 0 d0 50.7 60 OY) 


Fig. 3. — Caratteristiche del tubo E 80 T 
in particolari condizioni di lavoro. 


Con N di questi tubi in parallelo e con differenze di potenziale V, opportu- 
namente graduate tra le due placchette D e D’ di ognuno di essi, si ottiene 
un discriminatore differenziale a N canali ciascuno di larghezza 2AV. 


La buona risposta di questo 
tubo alle alte frequenze ci ha 
indotto a studiare, in primo 
luogo, la possibilità di utilizzare 
la E 80 T come discriminatore 
veloce. 

Elemento essenziale di un 
discriminatore di questo tipo è 
un formatore di impulsi rettan- 
golari, di ampiezza proporzio- 
nale all’ampiezza dell'impulso 
in esame, capace di funzionare 
correttamente con alte fre- 
quenze di ripetizione. 

Un primo modello di for- 
matore è disegnato in Fig. 4. 
Esso trasforma il breve impulso 
da analizzare in un impulso la 


10 - Il Nuovo Cimento. 
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Fig. 


. — Schema di principio del formatore 
di impulsi. 
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cui durata è controllata da una linea. I difetti di un tale sistema, (non per 
fetta linearità tra ingresso e uscita, frequenza non superiore a — 10 MHz) 
derivano essenzialmente dalla capacità presente tra la placca ed il catodo dei 
diodi. Ciò giustifica l’utilizzazione del pentodo E 80 L (Fig. 4) in sostituzione 
di un diodo che avrebbe in ogni caso una capacità pit elevata. 


Entrata 
Formatore | J L \ L 
di impulsi TOYO YX VOY TOO Y VOKO DOT “OTP OPO 
|proporzionale OV SV | OV] BV] 20v| 30v| 35V| 4ov 
K KI KI KI KR ZA 
da M Rosi : 
[A 4 A V, V, Vale V, kh [2 
2 3 4 5 6 7 8 9 + 
er OP ON OORG eee 
UE 
g2g1 
A /\ 
[a> 00000 90000 29000 00000 20000 09000 91000 09000 
L 
T.=200ns Ze 
-22V - -15V 


Fig. 5. — Schema di principio del discriminatore differenziale. 


Questo formatore deve anche fornire un impulso di accensione per le E 80 T 
durante la parte centrale dell'impulso in esame; nel nostro montaggio infatti 
le E 80 T sono normalmente interdette. 

La realizzazione di un discriminatore differenziale con E 80 T presenta 
essenzialmente una difficoltà: la struttura della E 80 T è infatti tale che quando 
il fascio di elettroni non cade sulla placca esso viene assorbito dalla placchetta 
a potenziale superiore. Questo effetto ‘provoca perturbazioni solo nei canali 
inferiori (*) a quello interessato, perchè solo in tali canali il fascio elettronico 
colpisce la placchetta D che riceve l’impulso in esame. 

Bastano però opportuni accoppiamenti induttivi tra le placchette D dei vari 
canali per evitare che tali perturbazioni giungano sulla placchetta della val. 
vola interessata in tempo per disturbarne il funzionamento. Ciò si realizza 
accoppiando le placchette D, ed anche le griglie di accensione g, degli N canali 
con linee di ritardo in maniera analoga a quanto si fa in un amplificatore di- 
stribuito. ‘ 

Se l’impulso su tali linee si propaga dall’ennesimo canale al primo, le per- 
turbazioni, che si generano solo nei canali inferiori a quello interessato, possono 
raggiungere questo solo con un certo ritardo rispetto al segnale di accensione. 


(*) Intendiamo per canali inferiori quelli destinati ad accettare gli impulsi di am- 
piezza inferiore all’ampiezza dell'impulso in esame. 
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Tale ritardo è evidentemente uguale al doppio del tempo di transito tra stadio 
e stadio delle linee di accoppiamento. Se quindi l’impulso di accensione ha du- 
rata inferiore a questo ritardo, le perturbazioni non possono disturbare il fun- 
zionamento del discriminatore. 

Tali perturbazioni potranno però interferire con l'impulso successivo a 
quello che le ha generate; la durata di un tale stato di perturbazione sulla linea 
è però non superiore al doppio del tempo totale di transito della linea. 

Con tale sistema si possono realizzare discriminatori con tempo morto non 
superiore a 0.5 us. Sono in fase di studio discriminatori molto più veloci, an- 
cora basati sull'impiego di E 80 T. 


Ringraziamo i dottori A. EGIpI, P. H. STOKER e M. ToLLER per la loro 
collaborazione. 


SUMMARY 


A new type of pulse-height analyzer using a beam deflection tube is described. 
The use of these tubes allows such a circuit to work at high repetition rate, without 
structural complexities. 
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Ricerca sulla esistenza di particelle 
di massa intorno a 550m, nella radiazione cosmica (*). 


M. ConveRrsI, G. M. DE MUNARI, A. EGIpI, E. FIORINI, S. RATTI, 
C. RUBBIA, C. Succi e G. TORELLI 


Istituti di Fisica delle Università di Milano e di Pisa 
Istituto Nazionale di Fisica Nucleare - Sezioni di Milano e di Pisa 


(ricevuto il 29 Dicembre 1958) 


Riassunto. — Sono riportati i risultati di una esperienza condotta allo 
scopo di controllare l’esistenza (ed eventualmente studiare le proprietà) 
delle particelle di massa pari a circa 550 m, recentemente segnalate da 
A. I. ALIHANIAN et al. Il risultato negativo ottenuto è confrontato con 
quelli più recenti di vari altri gruppi. Si discute in particolare la possi- 
bilità che le particelle in questione siano simulate da interazioni nucleari 
di mesoni x o di protoni. 


1. — Al Laboratorio Raggi Cosmici del Sabbione (1), a 2550 m s. m., si 
sta concludendo una ricerca delle particelle di massa intorno a 550 m, (*) 
segnalate da A. I. ALIHANIAN et al. (23). 


(*) Presentato al XLIV Congresso della Società Italiana di Fisica, Palermo 6-11 
Novembre 1958. 

(**) Indicata nel seguito con la dicitura « particelle 550 ». 

(1) E. FIORINI, C. GrorI, C. Succi e C. MorIGGIA: Suppl. Nuovo Cimento, 7, 129 
(1958). 

(2) A. I. ALIHANIAN et al. — Intern. Conference on High Energy Physics (Mosca, 
Maggio 1956). 

(8) A. I. ALIHANIAN, N. V. SHOSTAKOVICH, A. T. DADAIAN, N. V. FEDOROV e 
B. N. DERIAGIN: Zu. Eksp. Teor. Fyz., 31, 955 (1956); Soviet Phys. J.E.T.P., 4, 817 
(1957) e bibliografia ivi contenuta. 
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In una precedente comunicazione (*) sono stati descritti i due dispositivi 
sperimentali con cui furono eseguite due successive serie di misure (Fig. 1). 


Supporto 
di legno 


Fondo camera C.I 


__ Scintillatore Sc 


i 


Tetto 


Fondo 
camera CR. 


Scintllatore in anticoincidenza Sq 


Fig. 


=" camera CR. 


- 
—iC R\ 
zi 


LE 


Nella prima serie di misure sono stati raccolti ~ 2100 fotogrammi in 
~ 800 ore di effettiva misura; nella seconda serie ~ 3200 in — 900 ore di mi- 


sura. 


(4) M. ConveRrsi, G. M. DE Mounart, A. Earp, E. FIoRINI, S. RATTI, C. RUBBIA, 
C. Succi e G. TorELLI: Il dispositivo sperimentale per la ricerca di particelle di 550 me. 
Comunicazione al XLIV Congresso della Società Italiana di Fisica, Palermo, Novem- 


bre 1958. 
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Nelle nostre misure la massa delle particelle che si arrestano nella camera 
a lastre C R (Fig. 1) viene ricavata nella maggior parte dei casi dall’ampiezza 
dell’impulso fornito dallo scintillatore proporzionale S, e della misura del per- 
corso residuo. Quando l’impulso dello scintillatore S, non è utilizzabile (*) o 
quando comunque il risultato della misura sopra citata presenta qualche dubbio, 


dall’esame dei fotogrammi è sempre stato possibile identificare la particella, 
valutandone la massa mediante uno o più dei seguenti metodi: 


a) stima della ionizzazione in un punto della traccia e misura del per- 
corso residuo; 

b) valutazione della variazione di ionizzazione lungo la traccia nella 
camera C R; 

ce) misura degli angoli di « scattering » nelle lastre e del percorso residuo(*). 


Sono stati esclusi dall’analisi tutti i fotogrammi in cui la particella in esame 
termina il proprio percorso in una delle prime tre lastre della camera C R perchè, 
in questi casi, non si hanno sufficienti indicazioni per una sicura individuazione. 

L'applicazione di questi criteri si rende necessaria in particolare per rico- 
noscere quei protoni che, dopo aver attraversato lo scintillatore proporzionale 
S., subiscono una interazione nucleare prima di essere frenati in una delle 
lastre di C R. Poichè nella nostra esperienza la massa di una ‘particella viene 
determinata con misure di perdita di energia e di percorso residuo, occorre porre 
particolare attenzione ad eventuali interazioni nucleari dei protoni nelle lastre 
e nel tetto della camera C R; per effetto di tali interazioni un protone si pre- 
senta con un percorso residuo inferiore a quello che gli competerebbe se si fer- 
masse perdendo energia esclusivamente per ionizzazione e può così apparire 
con una massa inferiore a quella vera. Le interazioni nucleari in questione sono 
essenzialmente: 


a) interazioni tra lo scintillatore $, e la camera C R con un solo secon- 
dario visibile in C R; 
b) interazioni in una lastra di © R senza emissione di secondari visibili; 
e) interazioni in una lastra di C R con emissione di un solo secondario 
visibile. . 
In Fig. 2 è mostrata una interazione del tipo c). 
Nella prima serie di misure il riconoscimento delle interazioni di tipo a) 
è facilitato dall'esame della traccia nella camera di Wilson superiore (Fig. 1). 
Nella seconda serie di misure l’impulso del contatore di Cerenkov 0, unitamente 


(*) Per il 30 % circa del tempo utile di misura, lo scintillatore S, ha funzionato in 
modo non pienamente soddisfacente. I dati raccolti in queste condizioni sono stati 
analizzati a parte e inclusi successivamente nei risultati finali perchè trovati consistenti 
con gli altri. 

(5) M. Annis, H. $. BRIDGE e S. OLBERT: Phys. Rev., 89, 1216 (1953). 
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all'analisi della traccia in C R permette il riconoscimento di tutti i protoni che 
subiscono i tre tipi di interazione elencati. 


2. — I risultati ottenuti sono riportati nella Tab. I. Nelle prime due righe 
della Tabella sono riportati i numeri di protoni e di mesoni leggeri (4 e x) che 
si fermano per sole perdite di energia per ionizzazione in una delle lastre CR. 


TABELLA I. 

Percorso residuo in acqua 

Eventi da 59.1 a 64.3 g/em? | da 91.5 a 96.7 g/cm? 

14 serie di misure 2% serie di misure 
| Protoni a fine percorso . . . . . . 197 251 
Mesoni nw e 7 a fine percorso, . . . 5 5 
Protoni simulanti particelle 550 . . 9 16 
«pu equivalenti » rea ~ 1400 ~ 3600 

Particelle 550 calcolate in base ai 

dati di A. I. ALIHANIAN et al. . = 7 > 18 
| Particelle 550 osservate . . . . . . 0 0 


Nella terza riga della Tabella è riportato il numero di protoni che intera- 
gendo senza emissione di secondari visibili o con emissione di un solo secondario 
visibile, verrebbero scambiati per particelle 550 se la misura della massa fosse 
ricavata soltanto dai valori della perdita di energia nello scintillatore S, e del 
percorso residuo in CR. 

Nella quarta riga della Tabella è infine indicato, con l’abbreviazione «4 
equivalenti », il numero calcolato dei mesoni y che, entrando nell’angolo solido 
dell’apparato, sono ‘arrestati -in una delle lastre della camera CR durante la 
misura. 

Per il rendimento assai elevato del contatore C, di anticoincidenza, solo 
una frazione assai piccola di tali mesoni è stata effettivamente registrata (se- 
conda riga della Tab. I). 

Il numero di «x equivalenti » è stato ricavato indipendentemente in due 
modi diversi: 


a) estrapolando gli spettri di frequenza dei mesoni y alle varie quote (91°). 


(5) Progress in Cosmic Ray Physics, vol. I (Amsterdam, 1952), p. 340. 

(?) C. E. Mizter, J. E. HENDERSON, D. $. PoTTER, J. Topp jr., W. M. SANDSTRÒM, 
G. R. Garrison, W. R. Davis and F. M. CHARBONNER: Phys. Rev., 93, 590 (1954). 

(8) M. N. KoëarrAN: Soviet Phys. J.E.T.P., 1, 128 (1956); 3, 350 (1957). 

(9) B. Rossi: Rev. Mod.-Phys., 20, 537 (1948); B. Rossi, M. SAND e R. SARD: 
Phys. Rev., 72, 120 (1948). 

(19) M. Conversi: Phys. Rev., 79, 749 (1950). 
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b) moltiplicando l'intensità dei protoni ottenuta sperimentalmente per il 
rapporto tra il numero dei mesoni y e dei protoni fine percorso, ricavato dai 
dati di C. E. MILLER et al. (7). 


Se le particelle 550 avessero la frequenza riportata da A. I. ALIHANIAN 
et al., di una su duecento mesoni y, avremmo dovuto osservarne complessiva- 
mente circa 25: nessuna delle particelle rivelate con la nostra apparecchiatura; 
è stata riconosciuta come tale. 


3. — Il problema dell’esistenza delle particelle 550 è stato affrontato da 
numerosi altri autori con esperienze eseguite sia alle macchine sia coi raggi 
cosmici. 

Le esperienze alle macchine sono state eseguite per cercare le particelle 
550 intese come prodotti di decadimento o di interazione. In particolare: 


a) G. HARRIS et al. (2) hanno osservato in emulsione nucleare il decadi- 
mento di 5000 mesoni K* prodotti al Cosmotrone: nessuno dei prodotti di 
decadimento aveva massa prossima a 550 m,. 


b) L. W. ALVAREZ et al. (12) hanno studiato con una camera a bolle a 
idrogeno in campo magnetico, la possibilità che particelle 550, con una quantità 
di moto di 450 MeV/e e una vita media maggiore di 1.6:107$ s, fossero prodotte 
dal bevatrone. La sezione d’urto per produzione di tali particelle è risultata 
inferiore di almeno 300 volte a quella di produzione dei mesoni K~ generati con 
la stessa quantità di moto. 


c) P. C. STEIN (#) ha studiato l’eventuale produzione di particelle 550 
in Berillio da parte dei raggi y di bremsstrahlung dell’elettronsincrotrone di 
Cornell (energia da 550 MeV a 1.1 GeV). Nell’ipotesi che particelle di massa 
circa 550 m, siano prodotte in coppie con una vita media di 150 us, la sezione 
d’urto è risultata inferiore a 4-10-33 em?/nucleone. 


Le ricerche compiute nella radiazione cosmica sono elencate nella Tab. IX. 
Per brevità la discussione verrà svolta solo sulle esperienze di J. W. K EUFFEL 
et al., A. I. ALIHANIAN et al. e N. Basu e M. S. SINHA nelle quali si hanno 
indicazioni più o meno favorevoli all’esistenza di particelle 550. 


(+) G. Harris, J. OREAR and $. TAayLoR: Nuovo Cimento, 6, 1232 (1957). 

(1?) L. W. ALVAREZ, H. BRADNER, P. FALK-VAIRANT, G. D. Gow, A. H. ROSEN- 
FELD, F. T. SoLMITz e R. D. Trrep: Proc. Seventh Annual Rochester Conference, VIII. 
16 (April 1957). 

(13) P. C. STEIN: Phys. Rev. Letters, 1, 21 (1958). 
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1) Nell’esperienza di J. W. KEUFFEL et al. (14) le particelle in esame dopo 
aver attraversato uno scintillatore proporzionale $,, sono assorbite in un grande 
sintillatore S, completamente circondato da un altro scintillatore in anticoin- 
cidenza $,; un contatore di Cerenkov, pure in anticoincidenza, provvede ad 
escludere buona parte degli eventi corrispondenti a mesoni leggeri frenati in 
Ss. Lo spettro delle masse ricavate dalle ampiezze degli impulsi degli scintil- 
latori S, ed S,, presenta un picco intorno alle 550 m,. Gli Autori, non essendo 
in possesso di una tecnica di osservazione visiva dei singoli eventi, non sono 
però in grado di stabilire se le particelle da essi rivelate come particelle 550 
siano o meno simulate da protoni che subiscano interazioni nucleari in S,: Si 
conclude comunque che le particelle 550, se esistessero, avrebbero una frequenza 
almeno 500 volte inferiore a quella dei mesoni y nello stesso intervallo di per- 
corso residuo. 


2) Nell’esperienza di A. I: ALTHANIAN et al., già descritta nella precedente 
comunicazione (4), la misura della massa delle particelle osservate veniva otte- 
nuta da misure di curvatura e di percorso residuo. A tale proposito era già 
stata messa in luce la possibilità (1°) che tali eventi potessero essere attribuiti 
a mesoni 7 che, interagendo nelle lastre della camera di Wilson, dessero un 
percorso residuo apparente inferiore al proprio. Per minimizzare tale possibi- 
lità gli autori hanno compiuto misure di scattering e di variazione di ionizza- 
zione sulle tracce delle particelle nella camera di Wilson. Ci sembra però che, 
dato l'elevato spessore delle lastre ed il basso numero di «gaps » disponibili, 
le condizioni sperimentali fossero troppo sfavorevoli per compiere misure suffi- 
cientemente attendibili. 


3) N. BASU e M. S. SINHA (1°) in una esperienza sulla intensità dei mesoni 
u al livello del mare, eseguita con una camera di Wilson contenente 5 lastre di 
rame da mezzo pollice di spessore, hanno recentemente osservato 5 particelle 
cui attribuiscono una massa di circa 500 masse elettroniche in base a misure 
di ionizzazione, variazione di ionizzazione e percorso residuo. A nostro parere 
le informazioni che si possono trarre dall’esame della traccia di una particella 
che attraversa quattro lastre di così elevato spessore non possono permettere 


(14) R. L. Catt and J. W. KEUFFEL: Bull. Amer. Phys. Soc., 6, 321 (1957); J. W. 
KrurreL, R. L. Catt, W. H. SANDMANN e M. O. Larson: Phys. Rev. Letters, 1, 203. 
(1958); J. W. Keurrez, R. L. CALL e W. H. SANDMANN: Bull. Amer. Phys. Soc., 8, 
162 (1958). , 

(5) J. DAUDIN: Journ. Phys. et Rad., 14, 419 (1953); Compt. Rend. du Congr. de 
Bagnères de Bigorre, p. 244 (1953). 

(16) N. Basu e M. S. SinHa: Ind. Journ. Phys., 86, 259 (1958) e bibliografia ivi 
contenuta. 
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fi 
Autori È Quota Tecnica 
(in m s.m.) 
i hea Gs Spettrometro magnetico | 
A. I. ALIHANIAN et al. (2%) (Mosca) 3200 | camera a si 
J. OREAR et al. (17) (Columbia) 3 300 Emulsioni nucleari 
N. DureA PRASAD e 0. P. SHARMA i wine Jsioni foun 
(Bombay) (18) 3400 Emulsioni nucleari 
Fazio (M.I.T.) (19) 0 Emulsioni nucleari 
J. W. KEvFFEL et al. (14) (Utah) 1400 Elettronica 
| N. Basu, M. $. SINHA (15) (Calcutta) 0 Camera a lastre 
AH 1, (20). (Princet E Camera con campo magn 
. HENDEL et al. (7°) (Princeton) 3250 | e camera ale 
| ee ey pine Camera con campo magn 
M. Inoxt et al. (21) (Yamanashi Univ.) 0 con Taio 


Presente lavoro (??) (Milano-Pisa) 2550 Camera a lastre ed elettro 


una valutazione attendibile della massa. Inoltre ci sembra molto difficile che 
gli AA. abbiano potuto, in queste condizioni, riconoscere tutti i protoni simu- 
lanti per interazione nucleare una massa inferiore alla propria. D'altra parte 
va osservato che la curva di distribuzione delle misure di massa ottenute da 
determinazioni di ionizzazione e di percorso residuo presenta una lunga coda 
verso le masse maggiori dovuta alle fluttuazioni alla Landau della ionizzazione; 
di conseguenza le particelle 550 osservate potrebbero essere dovute anche a 
fluttuazione della misura della massa di mesoni u 0 x. 


(7) J. OREAR, E. BIERMAN e S. RosENDORFF: Postdeadline to the 1958 New York 
Meeting of the Amer. Phys. Soc. (Jan. 1958). 

(5) N. DurGA Prasap e 0. P. SHARMA: Second U.N. Intern. Conference on Peaceful 
Uses of A.E. (Geneva, Sept. 1958). 

(19) J. STEINBERGER et al.: CERN Symposium, 1958, p. 153 (1958). 

(2°) A. HENDEL, P. PrRoué e G. T. Reynorps: Bull. Amer. Phys. Soc., 8, 25 (1958). 

(21) M. Inoxi, T. YASAKI, M. Macurpa and Y. Marsuxawa: Phys. Rev., 105, 
1872 (1957). 

(22) M. ConveERsI, E. FIORINI, C. GIORI, C. Ruppia, C. Succi e G. TORELLI: Proc. 
Intern. Conference on Mesons and Recently Discovered Particles, XIII-47 (Padua- 
Venice, Settembre 1957); M. Conversi, E. Fiorini, S. RarTI, C. Ruppia, C. Succi 
e G. TORELLI: Nuovo Cimento, 9, 140 (1958). 
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Particelle Particelle 
Metodo di misura de 1 550 550 
ao previste osservate 
rvatura e percorso residuo 2200 — 11 
ndard 700 3—4 0 
ndard 500 2-3 0 
indard 1250 6 0 
‘dita di energia parziale e totale 28 000 140 < 50 
riazione di ionizz. e percorso residuo ? a 5 
rvatura e percorso residuo 1000 5 0 
rvatura e variazione di ionizzazione 80 1 < 3 
‘dita di energia e perc. residuo 5.000 25 0 


In conclusione noi riteniamo che in gran parte, se non tutti, gli eventi inter- 
pretati da A. I. ALTHANTAN et al. e da N. BASU e M. S. SINHA come particelle 
550 siano pioni o protoni che, prima di terminare il loro percorso, abbiano 
subito delle interazioni nucleari non identificabili con le tecniche impiegate. 

D'altra parte i risultati della presente esperienza, in cui il numero di infor- 
mazioni disponibili è sufficiente a riconoscere senza ambiguità tutte le intera- 
zioni nucleari dei protoni incidenti, portano a concludere che, se le particelle 
550 esistono, la loro frequenza è inferiore di almeno un fattore 10 a quella ri- 
portata da A. I. ALIHANIAN et al. 


SUMMARY 


The results are reported of a cosmic ray experiment carried out at an altitude of 
2550 m a.s.l. to search for the 550-mass particles reported by ALIHANIAN et al. No 
event has been found out of 25 expected on the basis of the rate given by these 
authors. This negative result is compared with the results of other workers. The pos- 
sibility for pions and protons to appear as 550-mass particles when undergoing uniden- 
tified nuclear interactions is discussed. 
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IL NUOVO CIMENTO Vou: XIN. 2 , 16 Aprile 1959 


Su un contatore di Cerenkov a gas ad alto rendimento (*). | 


M. BENEVENTANO (+), E. DE AGOSTINO (**), C. A. GALTIERI (**), | 
B. RISPOLI (**) (**) e A.-SERRA (**) 
(**) Istituto di Fisica dell’ Università - Roma 


(*"*) Istituto Nazionale di Fisica Nucleare - Sezione di Roma 
(*»*) Centro Coordinamento Elettronica del O.N.R.N. - Roma 


(ricevuto il 12 Gennaio 1959) 


Riassunto. — Si presentano i dati preliminari riguardanti il funziona- 
mento di un contatore di Cerenkov a gas. Le misure del rendimento del 
contatore, mediate sullo spettro dei mesoni p. dei raggi cosmici con quan- 
tità di moto maggiore di 1.75 GeV/c, in funzione della pressione del gas, 
hanno dato il valore di (96.7+1.1)% alla pressione di 10.23 atm N.P.T. 


1. — Introduzione. 


Uno dei problemi che si incontrano nelle esperienze eseguite all’elettrosin- 
crotone è la discriminazione fra pioni di momento di circa 1000 MeV/c, e quindi 
BP =0.990, ed elettroni di pari momento e B & 1. 

Allo scopo si è studiata la possibilità di realizzare un Cerenkov a gas con 
efficienza molto prossima all’unità (in modo da poter essere usato eventual- 
mente in anticoincidenza con altri rivelatori), pur mantenendo le dimensioni 
entro limiti ragionevoli (lunghezza circa 1 metro), e la pressione del gas (CO,) 
non molto elevata (< 11 atm), onde poter fare le pareti del contenitore suffi- 
cientemente sottili. i 

Perchè un contatore di Cerenkov a gas abbia le caratteristiche prima elen- 
cate, come è stato messo in evidenza da vari autori (12), è necessario disporre 


(*) Presentato al XLIV Congresso della Società Italiana di Fisica, Palermo, 6-11 
Novembre 1958. 

(1) A. AscoLI BALZANELLI e R. ASCOLI: Nuovo Cimento, 6, 1392 (1957). 

(2) R. J. Hanson e D. C. Moore: Nuovo Cimento, 4, 1558 (1956). 
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«di una catena elettronica capace di rivelare anche gli impulsi dovuti ad un 


singolo fotoelettrone emesso dal catodo del fotomoltiplicatore che rivela la luce 
della radiazione Cerenkov. 


‘2. — Deserizione del dispositivo sperimentale. 


Il disegno schematico del contatore di Cerenkov con le sue dimensioni prin- 
cipali è riportato in Fig. 1. 


| 900 mm 


ti 


Th) 


Pa | 


di 


SI 


PS 


Uno specchio sferico alluminato concentra la luce della radiazione Cerenkov 
sul catodo del fotomoltiplicatore montato sulla finestra A. La concentrazione 
della luce è favorita da una guida di plexiglas di 5 cm di spessore particolar- 
mente studiata. 

Questa disposizione geometrica si è resa necessaria per evitare gli inconve- 
nienti che si avrebbero se le particelle da rivelare attraversassero il catodo del 
fotomoltiplicatore. 

Per la misura dell’efficienza del contatore di Cerenkov, mediante i mesont y 
dei raggi cosmici, si è usato il dispositivo di Fig. 2. 

I contatori A e B sono due scintillatori plastici di forma circolare, di 17 cm 
di diametro ed 1 em di spessore. 
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La luce di scintillazione viene rivelata mediante due fotomoltiplicatori ad 
alto guadagno del tipo RCA 6810A. Gli impulsi da essi forniti hanno un’ampiezza 
sufficiente a comandare i circuiti di coincidenza che seguono. 
ZA Pb Il contatore C è uno scintillatore a toluolo di area 
=a 46x46 cm? realizzato nei-nostri laboratori. Quattro foto- 
moltiplicatori del tipo RCA 6342 rivelano la luce di scin- 

NS tillazione. 

I segnali forniti da tali fotomoltiplicatori vengono me- 
scolati a poi amplificati mediante un amplificatore rapido 
del tipo AR 62 realizzato dal Centro Coordinamento Elet- 
tronica del C.N.R.N. D è un fotomoltiplicatore Dumont 
6364 a grande area fotosensibile (5 in. di diametro) che rivela 
la luce della radiazione Cerenkov. I segnali da esso forniti 
vengono amplificati mediante due amplificatori rapidi del 
tipo AR 62 prima menzionati. 

L’assorbitore in piombo Pb, posto tra Be 0, e l’as- 


A 


C 3 sorbitore Pb’ impongono un limite inferiore all’impulso 
Fig. 2. (1.75 GeV/c) delle particelle rivelate dalla coincidenza 
tripla ABO. 


L’assorbitore Pb’ è stato studiato in modo tale da tagliare la componente 
molle dei raggi cosmici onde evitare un inutile affollamento nel contatore A. 

Per la misura dell’efficienza del contatore da noi studiato, la quale si riduce 
al conteggio delle coincidenze triple ABC e quadruple ABCD, ci siamo serviti 
della catena elettronica schematizzata a blocchi in Fig. 3. 


Amplificatore oincidenza 
rapida 


t=20ns 


Circuito 
di 
numeroz. 


La coincidenza ABC viene effettuata mediante un circuito di coincidenza 
tripla rapida (*) il cui tempo di risoluzione, in queste particolari condizioni di 
lavoro, è di 10 ns. 


(5) U. PELLEGRINI, B. RispoLi e A. SERRA: Nuovo Cimento, 9, 171 (1958). 


1070 


SU UN CONTATORE DI GERENKOV A GAS AD ALTO RENDIMENTO 159 


La coincidenza quadrupla ABCD è stata realizzata mediante un circuito di 
doppia con tempo di risoluzione di 20 ns. Un canale viene comandato dagli im- 
pulsi che provengono dal fotomoltiplicatore D e l’altro dai segnali forniti dal 
circuito di tripla ABO. 

Il circuito di coincidenza doppia lenta si è introdotto per evitare di contare 
degli impulsi di quadrupla senza che scattasse il numeratore di tripla come spesso 
accadeva (un simile inconveniente è stato riscontrato anche da HANSON (?) in 
un dispositivo sperimentale simile al nostro). 

Da quanto si è detto precedentemente si vede chiaramente che la parte più 
critica del nostro dispositivo sperimentale è la coincidenza doppia rapida che 
deve poter essere comandata anche dai segnali generati da un singolo elettrone 
emesso dal fotocatodo del fotomoltiplicatore D. (Un certo numero di impulsi 
da rivelare ha ampiezza comparabile con quella degli impulsi di fondo ed anche 
minore, poichè questi ultimi hanno una distribuzione (4) continua in ampiezza 
che va da 0 a & 8 fotoelettroni.) 

Per un corretto funzionamento bisogna quindi: 


a) amplificare sufficientemente gli impulsi provenienti dal fotomolti- 
plicatore D; 


b) ridurre entro limiti ragionevoli il numero delle casuali. 


(4) G. A. Morton e J. A. MicHELI: Nucleonics, 4, 1, 16 (1949). 
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Per lavorare in condizioni da soddisfare la condizione a) ed essere quindi, 
con buona approssimazione, in condizioni da rivelare anche un singolo foto- 
elettrone si è fatto uno studio preliminare del fotomoltiplicatore 6364. Lo spettro 
integrale degli impulsi di fondo (Fig. 4) ci ha permesso di conoscere oltre che la 
loro distribuzione in ampiezza anche il loro numero nelle varie condizioni di am- 
plificazione. I-risultati sperimentali son però poco precisi. 


4 E ia Tea 
x10 ne imp/s H LA ee Il Je | 
15 sla I FM 6364(Curva integrale) 
REI | IV.L.1400 V 
+ + TH + T | + 
| + + ia tot 
10 ot L 
| 4 pl: 
IE ee RE je) VI RO | 
| +++ IF + — + 
tes cal i aan 
SE a et eS bea a 
=” ets SI: LI 4 
| eee MELI Sal di a 
UE Ah) ee 
0 8 5 10 Men 
Fig. 5 


Ripetendo infatti più volte la stessa misura, ad intervalli di tempo, la 
curva integrale si spostava ‘parallelamente a se stessa pur mantenendo lo 
stesso andamento. Spesso inoltre si osserva, durante la stessa misura, un 
improvviso aumento del numero degli impulsi, non giustificabile con consi- 
derazioni statistiche. Questo effetto è stato anche osservato su un fotomol- 
tiplicatore del tipo RCA 6342. 

Attualmente non abbiamo dati sufficienti per formulare una ipotesi atten- 
dibile che spieghi i fenomeni prima messi in evidenza ma ci proponiamo di fare 
in seguito uno studio più accurato. 

Conosciuta la distribuzione integrale delle ampiezze degli impulsi, si è illu- 
minato il fotocatodo con una debole sorgente di luce e si è ricavata la distribu- 
zione spettrale degli impulsi nelle nuove condizioni di lavoro (Fig. 5). 

Si è potuto quindi facilmente ricavare lo spettro differenziale degli impulsi 
generati dalla sorgente di luce. 

La curva, (Fig. 6), presenta un massimo per un particolare valore V del 
taglio al discriminatore corrispondente alla ampiezza degli impulsi generati da 
un fotoelettrone. 
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_ Queste misure anche se non molto precise ci hanno permesso, come i risul- 
tati sperimentali confermano, di conoscere Vordine di grandezza dell’amplifi- 
cazione che bisognava appli- 

care agli impulsi provenienti ae 


TA T 
dal Cerenkov. ere ble] mea aia es 
È | JL | J 
Dallo spettro integrale ci |_|_|FM 6364(Curva differenziale) 
degli impulsi di fondo si vede + FE | | IL. 1400 V D 
onthe. 


che il loro numero è molto 
elevato (104~10°) il che 
potrebbe dare luogo ad un 
eccessivo numero di casuali. 
Per evitarlo possono adot- 
tarsi vari accorgimenti: 


1) ridurre al massimo 
il tempo di risoluzione del 
circuito di coincidenza; nel 
nostro caso non si può scen- 
dere al di sotto di 20 ns per- 
chè, date le notevoli dimen - 
sioni del fotocatodo, i fotoelettroni ché partono dai suoi vari punti non im- 
piegano lo stesso tempo per raggiungere l’anodo: le variazioni nel tempo di 
transito sono dell’ordine di 15-20 ns; 


2) ridurre al massimo il numero degli impulsi di fondo del fotomoltipli- 
catore. 


A questo fine oltre le normali precauzioni che un buon montaggio richiede 
si può operare come segue: 


a) raffreddare il fotocatodo (noi abbiamo evitato di farlo a causa delle 
complicazioni meccaniche che ciò richiede); 
b) evitare effetti di integrazione dei piccoli impulsi dovuti ai dinodi del 
fotomoltiplicatore (a questo scopo noi ci siamo serviti di amplificatori rapidi 
con tempo di salita 3.5 ns. Quest'ultimo accorgimento, a quanto riferito in (1), 
è di gran lunga il più efficace). 


3. — Risultati sperimentali. 


Nell'ipotesi che si sia in grado di rivelare impulsi dovuti alla radiazione 
Cerenkov che diano luogo anche ad un solo fotoelettrone, il rendimento del 
contatore Cerenkov è dato da: 


(1) F—1—Aexp[—N,(c)], 


11 - 71 Nuoro Cimento. 
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dove N,(c) è il numero medio di fotoelettroni emessi dal fotocatodo quando il 
contatore è attraversato da una particella con f & 1. 


0 1 = as 4 5 6 2 8 9 10 11 Press (atm) 
= = i = = = |Press. alla temp. di 0°C 
Za at Sei Tr = | t ] J i 
T TS A2 0, 
N CE =168% € 2407-10 *_ 
a He Sees 
10 


5555 === “bee 


E NNO COSC SEITEN 


Ete 


6.44 4 


10° 


ove: 


ns è Vefficienza del sistema ottico; 
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n; è Vefficienza di conversione fotoelettrica del fotocatodo del fotomoltiplica- 
tore; 
En =n—1 con n=indice di rifrazione del gas; 


I è la lunghezza del percorso delle particelle nel contatore. 


Riportando in un grafico semilogaritmico log (1 — F) in funzione di e, si 
ottengono delle rette di parametro 7;7:. Analogo risultato si ha nel caso che ven- 
gano rilevati due fotoelettroni emessi dal fotocatodo (Fig. 7). 

Dal confronto con i dati sperimentali si può quindi verificare se la condi- 
zione che un singolo fotoelettrone venga rilevato è soddisfatta. 

Nella Tab. I sono riportati i valori sperimentali F” delle efficienze alle varie 
pressioni e i valori F corretti per l’efficienza a contatore vuoto x=(5.79+0.44)%: 


TABELLA I. 


P, (atm) 


ni 7.12 | 8.36 | 9.35 | 10.23 


FS) 185.7 1.7 | 87.0-+1.5 | 91.9-1.1 | 96.141.0 | 96.9+0.9 | 96.941.0 | 


| | | | 
86.1+1.6 | 91.3412 | 95.8+1.1 | 96.7+1.0 | 96.7+1.1 | 


(%) | 84.8+1.8 


L’efficienza «, diversa da zero, si attribuisce a casuali ed essenzialmente a 
sciami estesi di raggi cosmici. 

Il pianerottolo che si nota alle alte pressioni è dovuto ad un non perfetto 
funzionamento dell’elettronica. Questo inconveniente è in fase di eliminazione. 
La bassa pendenza dei punti sperimentali indica, posto 7, = 0.10, un’efficienza 
per il sistema ottico 7,=~20%, giustificabile con il fatto che, a causa di 
alcuni incidenti di lavorazione, l’alluminatura dello specchio non è riuscita del 
tutto soddisfacente. 

Con la sostituzione dello specchio noi speriamo di migliorare notevolmente 
i risultati ottenuti. 


SUMMARY 


The performance of a CO, Cerenkov counter has been studied. The efficiency 
of the counter for u-mesons of cosmic rays with momentum greater than 1.75 GeV/c has 
been measured as a function of gas pressure. Preliminary results give for the effi- 
ciency the value (96.7 +1.1)% at the pressure of 10.23 atm N.P.T. 


12 - JI Nuovo Cimento. 
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LETTERE ALLA REDAZIONE 


(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lasciata 
dalla Direzione del periodico ai singoli autori) 


Note on Parity Violation in Strong Interactions. 


B. D'EspaGNAT and J. PRENTKI 
CERN - Geneva 


(ricevuto il 16 Marzo 1959) 


Recent experiments have led, as a 
very preliminary result, to the conjecture 
that parity might be violated in some 
of the strong interactions that involve 
strange particles (1), while it is con- 
served to a high degree of accuracy in 
nuclear physics phenomena (?). It has 
been pointed out that these facts seem 
difficult to reconcile with isotopic spin 
conservation: indeed if A is really a 
singlet in isospace the reported longitu- 
dinal polarization of the produced A par- 
ticles would be in contradiction with 
strict parity conservation in nuclear 
phenomena if the latter is a consequence 
of general symmetries alone. This was 
shown by Stapp (?). 

Under these circumstances it seems 
worthwhile to investigate whether a 
strong parity violation in strange par- 
ticle productions or interactions could be 
reconciled with some very weak and as 
yet unobserved parity violation in nu- 
clear phenomena. Such a quantitative 
question is of course difficult to answer: 
one may for instance very well believe 


(3) H. BLUMENFELD, W. CHINOWSKY and 
L. M. LEDERMANN: Nuovo Cimento, 8, 296 (1958). 

(3) See e.g. D. H. WILKINSON: Phys. Rev., 
109, 1603, 1610, 1614 (1958). 

(*) H. STapPP: to be published. 


that such a state of affairs results only 
from dynamical particularities (large 
energy denominators, velocity dependent 
parity mixtures, etc.). Calculations along 
these lines have not yet been done so 
that such a possibility still remains: 
Alternatively one may also believe that 
effects of this sort, although they prob- 
ably are present, are not sufficient to 
account for the surprising difference of 
behaviour of the two kinds of pheno- 
mena. Some general principles should 
then be sought by means of which 
parity violation should be preserved 
in, say, A production and quenched in 
nuclear phenomena. The purpose of this 
note is to sketch very briefly a possible 
way of approach in this direction. 

The general idea is first to split the 
strong interaction hamiltonian H,,, into 
two parts of different strengths, as has 
often been proposed in other connections: 


Ay = H,+ H,, 


where the very strong part H, could 
for instance be the 7 interactions and 
where the weaker part, H, could be 
the K interactions (*). It is assumed 


(*) In our example H,, contains only 
Yukawa type interactions. 
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that the coupling constants in H, are 
weak enough so that a perturbation 
expansion in H, (not in H,) makes some 
sense. H, is then himself split into 
two parts 


H, = gh+g'h', 


where À and h' could for instance be 
the NK and EK interactions. 

Let us next assume that H,, h and h’ 
can be chosen in such a way that there 
exist two distinct assignments of parities 
to the set of existing particle fields: 
one which makes H,+g% parity con- 
serving but not h', and one which makes 
H,+g'h' parity conserving but not h. 
Under these circumstances it is easy 
to choose H,, h and h’ in such a way 
that all matrix elements for nuclear 
reactions (*): 

a) do not contain any odd power 
oies 

b) are such that the terms of 
second order in H, do not involve the 
bilinear product hh’. 

For instance, conditions a) and b) 
are automatically satisfied if H,, h and h’ 
are the x, NK and EK interactions 
respectively. Then the parity violating 
terms in any nuclear matrix element 
are of order 4 or higher in H,: this we 
may assume will quench the parity 
mixture parameter in nuclear pheno- 
mena. 

Now we have still to choose our 
parity assignments in such a way that 
for A production no similar quenching 
should appear. For that purpose. 

i) we require that H,, h and h’ 
should all be invariant under the con- 
ventional charge symmetry operation 
(we do not require charge invariance). 
This condition is enough to make the 
self-mass matrix diagonal in A and XY, 
and guarantees therefore that the physi- 


(*) The term nuclear reaction will mean 
in the following all processes where the par- 
ticles in the initial and final state have all 
strangeness zero. 
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cally detectable particles are the quanta 
of the A and X, fields and not of any 
linear combinations of A and X, (**); 

ii) we choose, however, the assign- 
ments of parities to the various particle 
fields in such a way that the set of parity 
assignments which makes H,+ gh parity 
conserving involves in particular opposite 
parities for 


wi er 
Tn 
and 
pet 
V2 


This in no way contradicts any of 
our previous assumptions. It is true 
that, were h’ absent, condition ii) would 
imply that the self mass matrix is dia- 
gonal in Y and Z (this in turn would 
mean that Y and Z are the physically 
observable particles and together with i) 
would imply that A, 4), Y and Z have 
the same mass), but as long as g' 
is appreciable (g' > e) this is not the 
case. Thus nothing in condition ii) 
prevents A and %, from having different 
masses. 

h being the hamiltonian part which 
creates A+K from an initial system 
made of « ordinary » particles, the matrix 
element for A production, to lowest 
order in H,, involves h linearly (and H, 
to all powers). Because of the mixed 
parity character of the A this matrix 
element will in general give rise to an 
appreciable amount of parity mixture 


and, as is easily checked, to an 
appreciable longitudinal polarization 
of the A. 


Detailed hamiltonians that obey con- 
ditions i) and ii) can easily be written 
but considering the preliminary character 
of the experimental material this, up 


(**) As A and £ are transformed in A and 
— Y, under.charge symmetry, no transitions 
from A to ® or from to A states are pos- 
sible through H;nt- 


166 


to now, seems hardly worthwhile. A ge- 
neral and unpleasant feature of all these 
hamiltonians is that they do not conserve 
ordinary isospin: this is due to the fact 
that X, has a mixed parity while 2° 
and =~, being the only particles of their 
kind (specified Q, N and S), must 
necessarily be ascribed a definite parity 


B. D'ESPAGNAT and J. PRENTKI 


assignment. An approximate conserva- 
tion of isospin (in H,) can, however, 
always be introduced which might be 
sufficient for nuclear and pion physics 
phenomena. In such schemes the approx- 
imate equalities of the 2*, 2° masses 
would thus probably raise the most 
serious problems. 
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Sull’integrazione delle equazioni della magnetoidrostatica 
per configurazioni a simmetria assiale. 


E. CANOBBIO e R. CROCI 


Istituto Nazionale di Fisica Nucleare - Sezione di Milano 


(ricevuto il 20 Marzo 1959) 


Si espone brevemente un metodo per ricondurre l’integrazione delle equazioni 
della magnetoidrostatica, per un cilindro e un toro circolari, a quella di un’equa- 
zione integrodifferenziale non lineare. Per alcune condizioni al contorno richieste 
dal problema si possono assumere i valori sperimentali: in questo modo si tiene conto, 
nello studio dell’equilibrio delle forti scariche nei gas, di fenomeni, come l’inver- 
sione della componente assiale del campo magnetico, notati in alcune zone dei tubi. 

Le equazioni così ottenute si possono integrare numericamente, studiando anche 
la stabilità delle soluzioni rispetto a piccole variazioni delle condizioni al contorno; 
evidentemente questo è un aspetto, in pratica notevolmente interessante, del pro- 
blema più generale della stabilità magnetoidrodinamica. 


1. — In un sistema di unità di misura gaussiano razionalizzato, con c=1, le equa- 
zioni della magnetoidrostatica si scrivono: 


1 { grad p= (rot B)\B, 
(1) | divB = 0. 


In un sistema di coordinate ortogonali (x, x?, a), i coefficienti della cui metrica 
sono indipendenti da «, le soluzioni indipendenti da 3 soddisfano il sistema di 
equazioni: 


B B, OB. 
ey aap (E 0 ) | a3B; Sen 0, 


da! Ox? dal dal 
n Di À (ee ma) + SB, oe = 0, 

ang, hi + a? E NOI 

ait Ot gn faye vel +a®2B, a j pee SE 
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dove 


Va Van * dos “Ogg - 


Il sistema (2) è quasi-lineare di tipo composito; infatti le radici dell'equazione 
caratteristica sono: Sa 


a Bs = ahs Sa 
%=%=—y? 6,= Vall-q??, = —iV a -a®?. 
al! B, 
I sistemi di coordinate che consideriamo esplicitamente sono: 


a) le coordinate cilindriche; 


b) un sistema di coordinate (0, 0, g), legate a quelle cartesiane (x, y, 2) dalle 
relazioni 


| æ— À c089, 
{ y=Asng, A= R —0 cos 6, 
| 2 == 0 Sune! 


. 


valido nella regione toroidale 0<@<r, 0<0<27, 0< p<2a, r< R, dove R è il 
raggio maggiore del toro. I coefficienti della metrica sono: a!!=1, a??=1/9?, a®8=1/A?. 

In entrambi i sistemi considerati a!! è uguale ad 1. Moltiplicando allora la I (2) 
per B,, la II (2) per a?2B,, sommando e tenendo conto della III (2), si ottiene la 
equazione 


Op ap op > 


®) By oh + Bei = 0, 


che è formalmente identica alla III (2). 


2. — Sia «?=s(x', rt’) la famiglia di curve con parametro 7, definita dall’equa- 
zione differenziale 


Le formule di passaggio dal sistema (xt, x?) al sistema di coordinate (x’, 7’) e vice- 
versa, siano 

fa=a, fæ'= a, 

Il MRS (Te Il TS AO 


Nel sistema di coordinate (x, rt) le equazioni (3) e III (2) si scrivono: 


cioè p'= p'(t), Bi= Bit). 
Nel seguito apparirà come sia possibile determinare p'(r) e B;(t) quando siano 
assegnate le condizioni al contorno. 
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x 


3. — Per il seguito è utile introdurre il potenziale vettore A; nell'ipotesi che 
non dipenda da x, B, e B, sono dati da: 


1 04, 1 04, 


2 


À Va da’ DE Na. oui 


Derivando la I(2) rispetto a x?, la II(1) rispetto a x! e sottraendo, si ottiene: 


DO. 0B, 0 (OB, OB, 0B, OB,\ (0B, OB, 
4) B -| + a?B, > SERI 
AUX dx! | dx? = Ron (= | 1 ( dna È dx a dr | 
dai “da dBs 
È Boy ee et: 
GET KE OT 


4. — Conviene, d’ora innanzi, considerare separatamente le coordinate cilindriche 
e quelle toroidali. Tenendo conto della IV (2), in coordinate cilindriche la (4) 
si scrive: 


0 (0B, O0B,\ a?*B, 0 (dB, dB, 
do\ 0 do B, 00\00 do 


0 CEES, 0e Obs fa 
- (È (In va) + o a0 (In 4 a) (5; -3) 05 


Nel sistema di riferimento (x’, 7), ponendo (0B,/00 — 0B,/00) = 2: 


A Bs, 
—- ] en, 
da + A (In Va) 
cioè 
ex CZ ORA CUT) 
(5) ana; oppure a4, = 2, 
V 


dove A, è l'operatore di Laplace in due dimensioni, nel riferimento (x', 7). 
Per determinare c(t) è necessario assegnare t=S(o0, 9) per 0, 0 0, costante. 
L'equazione 


(6) —= 


h ower? n 7 x + z e, x ! 
è un’identità quando a 0 si sostituisca s(z’,t) e a 0, wv’. L’equazione (5) da 4;, 
e quindi Bj e Bj, in funzione di a’, +, s(x', t); sostituendone l’espressione in (6), 
sì ottiene per s(x’, 7) un’equazione cui imponiamo le condizioni al contorno: 
(7) (TE NCA Ch Qe = s(e',t(4)), 


dove (a) e t;(2') si possono dare ad arbitrio. Per esse, invertendo la (7) si ha: 


To(e) = S(9, 0), t1(0) = S(@, 2x). 
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Quando si conosca 4,4, per 0—0, è quindi possibile assegnare c(7); sia: 


Na) LOGIN 
“al $3 40. 


[AgAglo=0 = 


i 


Se To(e) = 0, 


È SN = Ble), 
o” Q 
da cui: c(t) = t°p(7). 
Se all’equazione 
I Ni et) 
data (21)? , 


imponiamo una condizione corrispondente nel sistema (+'7) alla: 


[Aslour a x(0) 0 


si ottiene formalmente: 


ne en aay Bln) à 
As (x ‘hance rs LO (Gr',, © E599) E — 8(& n) dé dy — 
Ero AIS di je s/0E\] ds 
-Z| EC mf gelo) — 5 Kw! vlE n) Cala MILL 


dove K è la fuzione di Green che soddisfa condizioni di Dirichlet omogenee per 
ld 
x =: 


2 


K(a', t|é, 9) = —21n (5) +2In (i) i 
fi P 


con 


1 
r= —V (w')? + E2— 2a’ E cos (s(a’, t) — 8(E, n)), 


Il 
Vo = Vw)? 62 + r? — 2u'E cos (s(x', t) — s(&, n)) - 
Ss 
L'equazione integrodifferenziale finale è la seguente: 


(8) 


0A, Os Os 1 (0A,0s ,,,04, ds 
Or da! OT (a)? \ Oa! Or Cori age 


By(r) si determina come e(r): se [Bs]o-o= y(0), si ha: 
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La condizione che sulla superficie plasma-vuoto la variazione di p+(B)?/2 sia nulla, 
può venire soddisfatta da una p opportuna, quali che siano B,, B,, B;, per o=7; 
per assegnare p(t) è però necessario aver integrato la (8), ottenendo s(x’, 7) e, quindi, 
S(o, 0). Inoltre si ha [B,|,_,=(d«(0)/d0) e, con Vipotesi che esternamente al plasma 
sia B,=0, «(0)=cost=[Ag],_,, è possibile assegnare f(g) se si assume [0B,/00]o_0=0. 
In tale caso si ha infatti: 


5. — Nel caso delle coordinate toroidali conviene sostituire all’equazione: 


dda; 
do B, 
il sistema di equazioni: 
| do 
— = a??B, 
| dt 
d 
| a si "i £ 
dt 


il cui integrale generale sia: 


Riuscendo ad esprimere B, e B, come funzioni di p(t, 7, 0) e v(t, 7, 0) il sistema: 


| = = B,(u(t, T; 0), v(t, T, 0), t, T, 0) , 


#1 00 
| de a? Ba (ut, v, 0), v(t, t, 0), t, t, 0) , 

con opportune condizioni al contorno definisce due funzioni 9=,(E, t, 0), 0—v(t, t, 0); 

per questo assegniamo il valore di y su una superficie cilindrica con generatrici 

parallele all’asse { e direttrice o = o(7); nel sistema di coordinate (t, 7), o è un 


parametro, e quindi si può porre o = o(t), così che il sistema (9) si riduce alla 
equazione: 


Dot, t)= a? B(u'(t,t), v'(t, 7), t,t), 


dove 
PAU T) Fr u(t, T, o(t)) e v(t, T) = v(t, qT; o(7)) . 


Nel sistema di riferimento che cosi abbiamo determinato, l’equazione (4) si scrive: 


CFA 10 = Oe Eira ee esl 
= +25 (In Va) = TR A te 
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dove 
On’ Ov! du Ov’ 
Cpuet ) Or to 


Integrando si ha: 


1 Va On, OF OX ep 
IR == pay i = 
ui ( i a 0 aa Ja +0) 


x 


Da qui innanzi il problema è analogo a quello del cilindro; eccetto per il fatto 
che l’equazione che si considera ricercando le soluzioni indipendenti dalla terza 
coordinata dell’equazione 4,2=f(0, 0) in coordinate toroidali, non è più un’equazione 
di Poisson. Non si è riusciti a calcolarne la funzione di Green; pensiamo che l’inte- 
grazione numerica si possa egualmente effettuare considerando il sistema di due 
equazioni differenziali. 

Per quanto riquarda il campo di integrazione, esso è dato dalle curve 


[ 4 = tu), 


pera. 


che devono soddisfare le condizioni 


| tO) ei, = ta ose 


Ti(W)u=u; = Tita(U)u=uy > 


(à = 1, 2, 3,4; per i= 4, à + 1— l; w è un parametro). 
Esse si determinano con le relazioni 


[ 0= (7), 
loin 


e imponendo, come condizioni al contorno per la funzione v'(t, 7), che sia 


Ringraziamo il prof. P. CALDIROLA per il suo costante interessamento a questo 
lavoro. 
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Sull’integrazione delle equazioni linearizzate della magnetoidrostatica. 


E. CaNoBBIO e R. CROCI 


Istituto Nazionale di Fisica Nucleare - Sezione di Milano 


(ricevuto il 20 Marzo 1959) 


In (1) si è mostrato un metodo per ricondurre l'integrazione delle equazioni 
della magnetoidrostatica, nel caso di un cilindro e un toro circolari, a quella di 
un’equazione integrodifferenziale non lineare. 

Vogliamo ora mostrare come il sistema che si ottiene linearizzando il precedente 
può essere ricondotto a un’equazione differenziale lineare. 

Si pensa che l’integrazione di questa equazione, e quindi il confronto, che sarebbe 
assai interessante, con le soluzioni del sistema della magnetoidrostatica (1), sia pos- 
sibile solo numericamente. 


1. — Se Be p, sono soluzioni del sistema (1) 


a) | grad p= (rot B)\B, 


divB = 0, 
corrispondenti alle condizioni al contorno del problema, perchè B+b e p,+yp sod- 


disfino ancora le (1) (al primo ordine in b e p), occorre che le perturbazioni b e p 
soddisfino il sistema linearizzato 


J grad p= (rot b)/\B + (rot B)\b, 
(2) | di 
ivb Qs, 


Sulla superficie fluido-vuoto devono sempre essere verificate le condizioni: 


f p+B-b)=0, 
(bn) i 0%, 


dove: m è il versore normale alla superficie, e con <x» si indica la variazione della 
grandezza x attraverso tale superficie. 


(1) E. CANOBBIO e R. Croci: Nuovo Cimento, 12, 167 (1959). 
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2. — Nelle coordinate cilindriche, e in quelle « toroidali » introdotte in (1), il 
sistema (2) si scrive: 


do, (2 BA) me, (O on, e 

a+ (Gig) rende (Ste) he 
(3) B, Li + a?°B, = = = bi = a22b,, 

È ace = = = (In Va)-b, — A (In 4/a@)+a22b, . 


Esso è di tipo composito, le radici dell’equazione caratteristica essendo 


Ba?” — 
0,= O2 = ; o3= iV a, o,= —iVa**. 


Ponendo 


= 33 
p = Gk, — Cd Bsrs , 


Dy = Cia, 
Dy = Cr, 
Ds = C525, 


dove €,, cs, C3, € sono funzioni arbitrarie di «1, x?, il sistema (3) diviene: 


Cy ie os À — Cs DI (B,a**)z3 — B,a??c, n a?2 a + Bae, Di de 
+ B,a?? A = a? Ce = a) Cy — 08863 = 23 
: oe 
6, = di 2 + ae, © aa PS, i = — Ca DI (in ae at ey À in veel 


Moltiplicando la II (3) per a?2B,/B,, sommando alla I (3) e ponendo 


sus oc 5 OB OBS ‘ LL 0G, 200, 1/0B, OB,\ _ 
i Mn UC rs EME Ca? da B,\Oa, da) ” 


1086 


— = 


SULL’INTEGRAZIONE DELLE EQUAZIONI LINEARIZZATE ECC. 175 


si ha 
0% WB, da CL DE a?°B, 0a? 22 OC B 
ca e 7 By Oat) ga OP? 
1/02, a? B, 0% Cia OB? ac, OB 
B;\ 0a B, 0a B, ox SB; On? “a> 
(4) : 
Ont (Va %) To (a224/a 2) = 0, 
dz CEA oe 38 
alee die 0, It 23 + a22B, Le 


dx! dx? Fm om ° da? 


Se s(x', t) è Vintegrale generale di dx?2/dx1 = a?*B,/B,, dal riferimento (21, x?) si passi 
a quello definito da (1): 

Ein 

Li SL, 7), 


nel quale le (4) si scrivono: 


0, , da'83 a 
= Za » 
Ox’ > Oa! pan 
| Oz, i : 
1 te 
x 
| asa DR NARRA na, 
Ot oa! (ava a By (Ca vVa'es) + ar (a’?°¢o/d'@) = 0, 
ds 02, ds da sa, 1021 08 Of, 08 Ode 
Se Ea ee A 122Ble 
Oc on! da Or a a à 
, (08 0a’ Os da'83 CCL 
= : 1222p! 4 
ee (È dd Oa! di | sera 
dove 
sta , dc, 0s/da' ci OB,  % OB, 0s/0x'\-, 
a'= a'?? a??' By — ; p= : ; 
dt Osfot B; 0x1 Bi dr Os/0r 


F ale, 0B, 1. \ 5; 
L Bi or > 4 


Dalle prime tre si ha: 
; da'83 
[0 = + as) da! + 10, 
23 - fa y'22) de' + g(t) , 


A/a'a 220,2, [5 (4/a' cat) DE i a cz) dr + h(x’) . 
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Ponendo 
1 — I os d CAE Os © MINE 
e Nara, T= — (4/a' cars) —— — (Vm C424) , 
€ Ox! dt Or Oa’ 


l’ultima si serive: 


Os Os OSO SARO 23) 
ah Si 3 d 
ale h+ ale (rar) aa ae ae Bt) T dt 
0s 0 08 da'83 à USE 
Ye d UA dy pny LS ha PT ) (4 
To (fear) ce al er ay fasi | NE aa! de (IE [ral a" Ja. 
ds © da'83 lo a 0 
7 d (| a Nes a a) dB, (6424) + a'2%5B; — = (eh) + 
La GB ds 0 Ta a'2?0! B es ; de a! Bi ds 0 Tar\+ ds da'8 4 
€ — h—— — —- 
TO” i n° De da 2 an! a] | dn! 
08 das ds da'3 ds da'33 ; 
tr i Bley da! + di frau È eo (|(ve{ 2 a" Jar) 
Nel caso del cilindro si ha, in forma pit compatta: 
SITO Ae , 08 0s 7 ; Os à 
Ean Caio a frar+ (è fra) dr'— a'°?cB5 23 Da (fra)+ 
4 2 (fra) = mere 


[On Ok. [a?2B,b, +.a8B,b, + pla-,= 0. 


dvi 


+ UE: AP. 


Le condizioni al contorno cui devono soddisfare d, e p sono 


Quelle per le 2; sono quindi: 
(5) p kee [a BC rano =O 


Nell’equazione per 24, 2; appare solo come @/0x'; una volta ottenuta 2,, 2, è nota 
a meno di una f(t) additiva, con la quale si soddisfa la (5). 


* * * 


Ringraziamo il prof. P. CALDIROLA per il suo costante interessamento a questo 
lavoro. 
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Nuclear Scattering of High Energy Electrons. 


E. PREDAZZI 


Istituto di Fisica dell’ Università - Torino 
Istituto Nazionale di Fisica Nucleare - Sezione di Torino 


(ricevuto il 31 Marzo 1959) 


Extensive numerical calculations of the diffraction patterns relative to a high 
energy electron’s collision with a coulombian field of a particular form, have been 
performed, recently (!). From these numerical calculations we have learned that 
in order to know the differential cross sections at large angles, an exact evaluation 
of the single partial wave phase shifts is needed. In fact for 236 MeV electrons 
seattered by gold, the differential cross section beyond 90° requires the knowledge 
of the phase shifts with an error not larger than 1075 rad. 

This fact has a simple theoretical interpretation: at high energies, small mo- 
mentum transfers are the most frequent, in our case, for instance, the differential 
cross section at 90° becomes about 10-4 times smaller than at 30°. This means 
that at large angles conspicuous interference phenomena must occur among the 
various partial waves which form the total scattered amplitude. For this reason 
the phase shifts need an extremely accurate evaluation. These interferential phe- 
nomena come out because of the large number of nodes involved in the high / partial 
waves which contribute to the scattering. For this one needs to remember that 
the function /sin0 P,(cos 0) may be approximated, for large 1, by a zero order 
Bessel function: 


4/Sin 0 P,(cos 0) ~ VIJSIM4+ 50], 


whose zeroes become more and more dense for large 1. 

The theoretical interpretation of the measured cross sections raises the problem 
of giving the partial wave phase shifts by means of a compact analytical form. 

Recently REIGNIER (2) elaborated an approximated method of solving a Klein- 
Gordon equation based on the solution by iteration of an equivalent Riccati’s 
equation. He obtained, in this way, values for the scattering phase shifts (for some 
charge distributions in gold) with an approximation of 5-10-* rad. 

We present here a more direct and simple approach, for the same problem, 
making use of a suitable simplification of the well known W.K.B. method. Electrons 


(* D: G. RAVENHALL and D. R. YENNIE: The Physical Society, 70, 857 (1957). 
(2) J. REIGNIER: Nucl. Phys., 7, 526 (1958). 
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are still considered Klein-Gordon particles as for Ravenhall-Yennie and Reignier. 
This is, in fact, allowable at high energies, since geometrical optics is expected to 
become valid. 
We effected a numerical control by computing the differential cross section 
for 236 MeV electrons which are scattered by a shell charge distribution in gold. 
The Klein-Gordon equation reads: 


du(r) | +3) + $) 


x = LE — V@MP+ m4] u(r) = 0, 


(1) 


where V(r) is the usual electrostatic potential. Putting: 
U(r) = V(r)(2# — Vir)) , 


the well-known W.K.B. phase shifts are: 


1 1)2 po l 1)2 
(2) m= | Ve Dir) ti dr il pe LE ae, 
T LR 


where 7) and r, are the turning points of the respective functions under the sign 
of integration. 

Neglecting at first the difference between 7) and r, through a first order expansion, 
we obtain an approximated expression of (2): 


Tai o 
(3) n= — sl V (\ Z2+ na de, with ke=/(kr)?—(1+ 32. 


0 


This formula however is not sufficiently accurate, but a considerable improvement 
is reached if the expansion includes terms of second order only. 
Since the internal and external potential for our charge distribution are: 


Ze? 

Vie. (< R) 
Ze? 

ae (r> È) 


the phase shifts with KR >1+44 are given by the formula: 


_ Ze i A C+? 
(4) m= bali NBA MIRE ET 


i (+ =) LEA Es 2) 
E } (2kR)?| 3kB 


Ze? 
i| + kia (liete) 5 
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while those relative to angular momenta 1+4> KR are: 


1 NZ 1 oe A 
Peters EN 


Ze? 
(5) pe aa) 


Obviously an approximation of this kind can be made for an arbitrary charge 
distribution. We have limited ourselves to consider the particular case of a shell 
charge distribution in order to compare our numerical values with those given by 
RAVENHALL and YENNIE (!). The numerical results corresponding to kR=5.6 are 
listed in Table I. 


Tape I. — Phase shifts for scattering from the shell charge distribution for gold. 
l ne ™ l i | M | 
0 — 0.820 06 — 0.81847 6 — 1.08480 | — 1.08313 
1 — 0.83253 — 0.829 82 7 Oso II 6014 
2 — 0.855 21 — 0.852 61 8 — 1.23632 | — 1.23679 
3 — 0.88041 — 0.88095 9 —1.30020 | —1.30053 
d — 0.926 64 — 0.92470 10 — 1.357 72 | — 1.357 94 
5 — 1.000 54 — 0.997 82 11 AO ON | S EN- ATOS, 


n are the phase shifts of RAVENHALL and YENNIE (1), 7, are the phase shifts 
corresponding to our approximation (eqs. (4), (5)). 

From Table I we see that the values so obtained are good with an error not 
larger than 3-10-? rad. 

The differential cross section values for 236 MeV electrons scattered by gold 
can now be evaluated making use of the so-called «reduced series » introduced by 
RAVENHALL, YENNIE and WiILson (5) for the scattering amplitude, and our phase 
shifts. The results are given in Table II. 


TaBLe II. — Numerical results for differential cross section. 


x 


) 40° 50° 60° 70° 80° 90° 100° 


3) |3.21386 - 10-27/8.24867 - 10-28 | 3.772 37 -10-° 19.190 74 -10-29 4.32761 - 10-2” |4.20327- 10-*°| 4.76364: 10-99 


These values, compared with the curve reported by HorsTADTER (4), show that, 
for angles not larger than 70°--80°, the agreement is rather good, the relative error 
being less than 10%. 

It is not possible to use our approximation, which implies a phase shift error 
of the order of 10-2 rad, at angles larger than 70°--80°. We wish to observe, 
however that a better knowledge of the first two phase shifts is sufficient for a 


(3) D. G. RAVENHALL, D. R. YENNIE and D. G. WILSON: Phys. Rev, 95, 500 (1954). 
(4) R. HOFSTADTER: Annual Review of Nuclear Science, 7, 231 (1957). 


180 E. PREDAZZI 


good evaluation of the differential cross section also at 90°. For this one has, for 
instance, the asymptotical values of the phase shifts valid for the first angular 
momenta of Verde (5). The theoretical interpretation of the forward differential 
cross section (9< 90°) for which our approximation is valid serves already to diseri- 
minate among the different possible charge distributions. 

Besides, the absence of the spin interaction in all the above calculations puts 
in discussion the comparison between theoretical and experimental values for 
angles > 90° being, for these angles, the differential cross section sensitive to single 
phase shifts variations of the order of 10-4 rad. 


** * 


I wish to thank Professor M. VERDE for having suggested this problem and 
for many helpful and interesting discussions. 


(5) M. VERDE: Nuovo Cimento, 10, 560 (1958). 
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